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Abstract 


This  thesis  describes  an  investigation  into  the  emission  properties  of  Resonant  Cavity 
Light  Emitting  Diodes  (RCLED).  RCLEDs  are  a  new  photonic  device  that  promise  laser-like 
properties  of  narrow  output  frequency  spread  and  high  output  power,  without  the  usual 
laser  disadvantages  of  shorter  lifetime  and  increased  maintenance  and  lifetime  costs. 
Specifically,  this  thesis  examines  the  functional  dependence  between  emitted  wavelength 
and  the  RCLED  microcavity,  and  attempts  to  correlate  theory  with  experimental  data.  This 
thesis  introduces  a  classical  wave  interference  model  that  calculates  the  Spontaneous 
Emission  Enhancement  Factor  (S)  with  angular  dependence. 
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Design  and  Characterization  of  Resonant  Cavity  Light  Emitting  Diodes 

1.  Introduction 

1.1  Background 

Light-emitting  diodes  (LEDs)  have  long  been  a  trusty  workhorse  of  the  scientific  and 
industrial  communities.  Though  relegated  to  display  and  short  distance  communications 
applications,  they  have  nevertheless  performed  well  and  faithfully.  The  simplest  LEDs  are 
composed  of  semiconductor  material,  such  as  GaAs,  with  a  p-n  homojunction.  When  an 
applied  forward  bias  exceeds  the  energy  bandgap  of  the  material,  injected  electrons  and 
holes  recombine  at  the  p-n  junction  to  form  photons  resulting  in  emitted  light.  The 
wavelength  of  the  emitted  light  (Jig)  is  related  to  the  energy  bandgap  by  E  =  he  /  A0 ,  where 

E  is  the  energy  bandgap  in  electron-volts,  h  is  Planck’s  constant,  and  c  is  the  speed  of  light 
in  a  vacuum.  The  failings  of  LEDs  had  always  been  low  intensity,  short  coherence  time, 
and  omni-directional  output.  Thanks  to  recent  advances  in  the  epitaxial  growth  of 
compound  semiconductors,  LEDs  are  now  poised  to  enter  a  new  era.  Newer  LEDs  have  the 
potential  to  replace  incandescent  bulbs,  and  even  laser  diodes,  in  a  host  of  new  applications 
such  as:  automobile  brakelights,  flat-panel  high- definition  displays,  holographic 
illumination,  barcode  scanning,  and  as  light  sources  for  optical  interconnects,  fiber-optic 
communications  networks,  and  optical  data  storage  systems  [1;  2;  3].  The  useful  lifetime  of 
LEDs  (>  1  million  hours)  also  promises  low  replacement  costs  in  these  applications  [4], 

The  topic  of  this  thesis  is  a  new  architecture  that  promises  LEDs  with  comparable 
coherence  time,  wavelength  spread,  and  beam  characteristics  to  lasers.  This  architecture  is 
known  as  the  Resonant  Cavity  (or  Vertical  Cavity)  Light-Emitting  Diode  [5]. 
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The  Resonant  Cavity  Light-Emitting  Diode  (RCLED;  pronounced  rik-led)  has  many 
advantages  compared  to  conventional  LEDs.  Architecturally  it  is  very  similar  to  the 
Vertical  Cavity  Surface  Emitting  Laser  (VCSEL  pronounced  vik-sell)  in  that  it  consists  of  a 
quantum  well  active  region  in  a  Fabry- Perot  microcavity  of  width  d-q/ 1/2,  where  q  is  a 
positive  integer  (Figure  1).  Purcell  predicted  that  the  spontaneous  emission  from  a  medium 
within  a  microcavity  (mirror  separation  d  is  comparable  to  the  emission  wavelength)  would 
experience  a  change  in  lifetime  [6].  Schubert  et  al.  predicted  that  the  spontaneous  lifetime 
is  shortened  (enhanced)  for  on-resonance  optical  transitions,  but  is  lengthened  (inhibited) 
for  off-resonance  optical  transitions  [5].  The  RCLED  is  designed  so  that  the  spontaneous 
emission  lifetime  is  shortened  and  photons  are  emitted  at  an  increased  rate  with  smaller 
wavelength  spread  at  the  design  wavelength,  X0  [7;  8;  9].  Light  emission  occurs  through 

both  top  and  bottom  surfaces,  but  is  predominantly  through  the  mirror  with  the  lesser 
reflectance;  this  is  usually  designed  to  be  the  upper  mirror  ( Ripper  <  Riower).  Figure  1  shows 
the  differences  in  architecture  and  photon  emission  between  the  conventional  LED  and  the 
RCLED.  Figure  1  also  shows  the  angular  dependence  of  the  emitted  intensity.  The 
RCLED’s  emission  spectrum  changes  with  increasing  angle  from  the  surface  normal.  This 
is  not  the  case  for  the  conventional  LED.  Figure  2  shows  representative  emission  spectra 
vs.  relative  intensity  for  Figure  1.  The  emission  wavelength  and  full-width  at  half 
maximum  (FWHM)  of  a  RCLED  are  variable  and  are  controlled  by  the  cavity  width  and 
reflectances  of  the  upper  and  lower  mirrors. 

1.1.1  Distributed  Bragg  Reflectors.  A  DBR  consists  of  alternating  layers  of  high  and 
low  refractive  index  semiconductor  epitaxial  layers.  Each  DBR  interface  reflects  a  fraction 
of  the  incident  light.  For  the  DBR  to  be  strongly  reflective  the  light  must  add  up  in  phase, 
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and  therefore  each  layer  should  be  a  quarter  wavelength  (/i/4)1  thick.  Figure  3  shows 
the  alternating  high-low  (HL)  index  periods  of  the  DBR.  The  design  wavelength  (Xf)  of  the 

DBR  is  referred  to  as  the  Bragg  Wavelength.  The  DBR  Reflectance  is  a  function  of 
wavelength,  and  is  a  maximum  at  the  design  wavelength  at  zero  degrees  incident  angle. 
Maximum  reflectance  increases  with  increasing  number  of  HL  periods  and  also  increases 
for  larger  differences  between  the  refractive  indices  of  the  high  and  low  index  layers. 

1.1.2  Fabry-Perot  Cavity.  A  Fabry-Perot  cavity  is  formed  by  placing  two  mirrors  a 
known  distance,  d,  apart  and  ensuring  they  are  parallel  with  a  high  degree  of  accuracy  [10: 
159].  The  cavity  supports  longitudinal  modes  where  d-qX!  2 ,  where  q  is  a  positive  integer 
and  X  is  the  mode  wavelength.  The  design  wavelength  of  the  RCLED  must  be  a  mode  of  the 
cavity.  The  Fabry-Perot  cavity  of  the  RCLED  contains  a  photon  emitting  active  region,  for 
example  a  quantum  well  (QW).  The  quantum  well  must  be  designed  to  emit  photons  of  the 
proper  wavelength  to  couple  with  the  transmission  modes  of  the  Fabry-Perot  cavity.  The 
cavity  is  designed  so  that  transmission  modes  are  sufficiently  far  apart  in  wavelength  such 
that  the  medium  can  only  couple  (and  therefore  emit)  to  a  single  mode  within  the  QW 
spectral  gain.  Figure  4  shows  a  Fabry-Perot  cavity  with  incoming  (i),  reflected(r)  and 
transmitted©  electric  fields.  The  cavity  transmits  those  wavelengths  that  correspond  to  it’s 
longitudinal  modes. 


1  Note  that  X  =  X0/n,  where  n  =  refractive  index. 
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/  l  X 


Regular  LED 


/  l  X 


Resonant  Cavity  LED 


Note:  Length  of  arrow  denotes  relative  intensity  of  emitted  photons 


Figure  1.  Simplified  structure  of  LED  and  RCLED 
architectures  showing  photon  emission  patterns 


Figure  2.  Comparison  of  RCLED  and  conventional  LED  emission 
spectra 
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Incoming  Light 


Reflected  Light 


^(/4% 

^(/4nL 


Figure  3.  DBR  stack  containing  high-low  index  quarter  wave 
layers 


Ml  M2 

Figure  4.  Typical  Fabry-Perot  Cavity 


1.1.3  Quantum  Wells.  Quantum  well  (QW)  active  regions  lower  threshold  current  and 
simplify  wavelength  selection  in  laser  diodes  [11].  In  contrast,  LEDs  are  thresholdless  but 
can  also  use  QWs  for  more  efficient  photon  generation.  Figure  5  shows  the  conduction  and 
valence  bands  for  a  single  QW,  along  with  the  variables  that  define  the  quantized  energy 
levels  within  the  QW.  These  are  LQW:  the  width  of  the  well,  AEc:  the  conduction  band  offset, 


and  AE  :  the  valence  band  offset. 
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Figure  5.  Simplified  energy  band  diagram  of  a  quantum  well 
heterostructure  showing  allowed  discrete  energy  transitions 
between  electron-hole  levels 


1.1.4  Theory  of  Operation.  Figure  6  shows  a  schematic  of  an  example  fabricated 
infrared  (IR)  RCLED  [15:  124].  The  upper  and  lower  mirror  structures  are  composed  of 
AlAs/Al0 15Ga0  85As  distributed  Bragg  reflectors  (DBRs).  The  applied  voltage,  Vin,  causes  the 
carriers  to  migrate  toward,  and  recombine  in,  the  quantum  well.  The  Fabry-Perot  cavity  is 
resonant  for  recombination  wavelengths  that  correspond  to  cavity  transmission  modes.  The 
spontaneous  emission  is  enhanced  for  those  modes,  given  that  the  QW  layer  has  been 
located  at  an  anti-node  of  the  electric  field  in  the  cavity. 

Band  offsets  are  barriers  to  current  flow  and  increase  device  series  resistance.  Schemes 
to  lower  the  series  resistance  include  compositional  grading  of  the  V 4  layer  interfaces,  5- 
doping  the  a/4  heterolayers,  using  non-AlxGa1.xAs  materials  for  the  DBRs  such  as 
ZnSe/CaF2,  and  straining  the  DBR  by  mismatching  the  lattice  constants  of  the  X/4  layers. 
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[12;  13;  15].  This  thesis  deals  exclusively  with  compositionally  graded  AlxGal  xAs  DBR 
designs. 

To  sum  up,  designing  the  RCLED  requires  matching  the  allowed  QW  energy  transitions 
with  the  transition  modes  of  the  Fabry-Perot  cavity,  and  placing  the  QW  at  an  anti-node  of 
the  cavity  electric  field  [9;  15]. 


Ge/Au/Ni/Au  n  contact 


Figure  6  Schematic  diagram  of  RCLED  test  structure 


1.1.5  Performance  Characteristics.  The  performance  advantages  of  an  RCLED 
compared  to  conventional  LEDs  are  that  the  RCLED  combines  the  thresholdless  operation 
of  an  LED  with  laser-like  coherence-time  and  beam  properties.  A  conventional,  bulk 
homojunction,  LED’s  full-width  at  half-maximum  (FWHM)  may  be  approximated  by 
AX=1.45(lP2kBT),  where  kB  is  Boltzmann’s  constant  and  T  is  the  temperature  in  Kelvin  [28: 
584;  5].  Thus  an  LED  emitting  at  850nm  at  room  temperature  (T  =  300K)  will  have  a 
wavelength  spread  (AX)  of  -  28 nm.  In  contrast,  an  RCLED  emitting  at  the  same  wavelength 
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could  have  4 nm  <  Aa<  28 nm  depending  on  the  mirror  reflectances  (see  experimental 
results)  [5]. 


The  benefits  of  a  smaller  AX  may  be  seen  from  a  fiber-optic  communications  example 
[14].  The  pulse  broadening  due  to  chromatic  dispersion  (caused  by  AX)  is  given  by: 


At  ■ 


l^o >W  (sec) 

c  dXn  ) 


(1) 


where  l  is  the  fiber  length,  AX  is  the  source  spectral  width,  and  the  bracketed  term  is 
the  chromatic  dispersion,  which  is  constant  for  a  given  fiber  (~100 ps/(km  nm)  at  875 nm  for 
silica).  Thus  by  decreasing  AX  by  a  factor  of  7,  pulse  broadening  is  reduced,  and 
consequently  the  length  between  optical-repeaters  in  the  fiber  is  increased  by  a  factor  of  7. 

1.2  Problem  Statement 

This  thesis  investigates  the  spectral  gain  characteristics  of  RCLEDs.  This  includes  an 
investigation  of  the  functional  dependence  between  emitted-intensity,  wavelength, 
radiation  pattern,  and  device  structure.  Achieving  this  goal  required  designing,  fabricating, 
and  testing  RCLEDs. 

Designing  the  RCLED  required  first  choosing  a  desired  emission  wavelength,  X0  and 

then  “tuning”  the  quantum  wells,  Fabry-Perot  cavity  and  DBRs  to  that  wavelength.  Tuning 
the  QW  involves  selecting  the  composition  of  AlxGa1.xAs  for  the  barrier  and  well,  and  the 
well  width  necessary  for  electron-hole  recombinations  at  the  design  wavelength.  Tuning  the 
Fabry-Perot  cavity  to  the  design  wavelength  requires  making  it  an  integral  number  of  XJ2 
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thick  at  Tuning  the  DBR  requires  that  each  HL  period  be  XJ2  thick  at  A,0.  These 
thickness’  are  specified  only  at  X0  due  to  the  wavelength  dependence  of  the  refractive  index. 

The  devices  were  grown  and  fabricated  by  the  Center  for  High  Technology  Materials 
(CHTM)  at  The  University  of  New  Mexico  (UNM)  at  Albuquerque  NM. 

Testing  the  devices  involved  electro-luminescence  (EL)  and  photo-luminescence  (PL) 
excitation.  The  resulting  emission  from  the  RCLEDs  was  sampled  for  relative  intensity  and 
wavelength  vs.  increasing  angle  from  the  normal.  Reflectance  spectra  were  taken  from  each 
sample  to  verify  the  reflectances  calculated  during  the  design  phase  and  also  measure 
effects  of  mismatches  between  gain  spectrum  and  the  Fabry-Perot  resonance.  A  model  was 
developed  to  apply  classical  wave  interference  to  the  cavity. 

1.3  Scope 

The  scope  of  this  thesis  is  to  design  and  characterize  the  spectral  properties  of  infrared 
(IR)  (to  ~  850 nm)  AlxGa1.xAs  RCLEDs.  The  actual  research  proceeded  in  the  following  steps. 

a.  Design  RCLEDs  with  numerical  software  using  previously  proven  techniques 
[15;  10;  20]. 

b.  Estimate  the  spectral  characteristics  of  the  RCLED  with  particular  emphasis  on 
the  functional  dependence  between  emitted  intensity,  peak  wavelength,  and  the  radiation 
pattern  [8],  Calculate  the  expected  spectral  characteristics  of  the  RCLED  for: 

Structure  1:  no  DBRs  around  the  resonant  cavity; 

Structure  2:  two  DBRs,  with  reflectances  Rlower~90°/»  and  RUpper~800/°> 
below  and  above  the  resonant  cavity  [16;  17]. 
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c.  Construct  the  two  RCLEDs  mentioned  in  b. 


d.  Characterize  emission  radiation  of  each  RCLED  using  photo-luminescence  (PL), 
and  electro-luminescence  (EL)  [18;  15]. 

Structure  1 


Figure  7.  The  two  structures  examined 

1.4  Approach/Methodology 

The  following  theoretical  basis  was  used  for  modeling  the  RCLED.  First,  the  allowed 
energy  states  in  the  quantum  wells  were  used  to  calculate  a  gain  spectral  model,  in  order  to 
get  the  expected  EL  of  Structure  1  (resonant  cavity  only)  [16;  19].  This  model  gave 
information  about  the  gain  characteristics  and  peak  emission  wavelength  of  Structure  1. 
Secondly,  the  Macleod  method  was  used  to  generate  the  reflectance  values  for  both  a  DBR, 
and  a  Fabry-Perot  cavity  formed  by  two  DBRs  [10;  20].  This  second  model  included  the  non¬ 
normal  propagation  of  light  and  absorbing  quarter-wave  layers 

1.4.1  Design  Assumptions.  Three  key  assumptions  are  used  as  a  basis  to  begin  work 
on  this  proposal.  First,  that  the  A./4  thick  DBR  layers  can  be  built  within  allowed  tolerances. 
Second,  that  the  composition  of  the  quantum  wells  will  be  accurate  enough  for  calculation 
of  material  properties  using  parametric  equations.  These  first  two  assumptions  are  based 


Structure  2 
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on  previous  demonstrations  of  technical  excellence  in  this  field  by  The  University  of  New 
Mexico  at  Albuquerque.  Third,  the  complex  refractive  index  information  available  for 
A1  Ga,  As  is  considered  accurate,  and  no  further  exploration  of  that  work  will  be  required. 

The  different  refractive  indices  of  compound  materials  are  published  as  a  function  of  the 
fractional  composition  variable,  x  [21;22]. 

1.5  Sequence  of  Presentation 

Chapter  2  will  familiarize  the  reader  with  the  physics  behind  each  of  the  components  of 
the  RCLED.  That  is:  the  compound  semiconductor  materials,  complex  refractive  index, 
distributed  Bragg  reflectors,  Fabry- Perot  cavities,  and  quantum  wells.  Chapter  2  will 
conclude  by  putting  these  elements  together  and  forming  a  RCLED.  Chapter  3  will  present 
the  RCLED  designs.  This  will  include  complete  simulation  of  the  RCLED  characteristics  for 
the  designed  structures.  Chapter  4  contains  the  experimental  setups  that  were  used  and  all 
of  the  collected  data  for  each  of  the  RCLED  structures.  Chapter  5  presents  the 
experimental  results  and  compares  them  to  the  theoretical  design.  Chapter  6  contains  a 
summary  of  conclusions  and  suggestions  for  further  research. 
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2.  Design  Theory 


This  Chapter  contains: 

•  the  physics  of  the  individual  components  of  the  Resonant  Cavity  LED  (RCLED). 

•  design  criteria  for  combining  the  individual  components  into  an  RCLED  structure. 

•  a  classical  wave  interference  model  for  calculating  the  spectral  output  of  an  RCLED. 

The  basic  components  of  the  RCLED  are:  1.  the  distributed  Bragg  reflectors;  2.  the  Fabry- 
Perot  optical  cavity;  and  3.  the  quantum  well  active  region  embedded  in  the  cavity.  Each 
sub-section  discusses,  in  some  detail,  the  nuances  behind  decision  making  in  RCLED  design 
and  how  each  individual  piece  of  the  RCLED  must  be  designed  and  built  to  work  at  the 
design  wavelength.  The  critical  parameter  is  the  design  wavelength,  or  peak  emission 
wavelength. 

2.1  Distributed  Bragg  Reflectors  (DBRs) 

The  distributed  Bragg  reflector  consists  of  periodic,  alternating  high-low  (HL)  quarter 
wave  (X0/4)  layers  ofrefractive  index  material.  The  layers  are  constrained  to  have  an  optical 

thickness  of  a  quarter  of  the  design  wavelength: 

l HnH  =  lLnL  =  A0  /  4  (nm)  (2) 

•  Ih  is  the  physical  thickness  of  the  high  index  layer 

•  II  is  the  physical  thickness  of  the  low  index  layer 

•  rw  is  the  refractive  index  of  the  high  index  layer 

•  til  is  the  refractive  index  of  the  low  index  layer 

•  Ao  is  the  design  wavelength  of  the  DBR. 
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The  quarter  wave  (Ao/4)  criterion  for  layer  thickness  ensures  that  at  normal  incidence, 
and  for  the  design  wavelength,  each  reflection  from  subsequent  layers  interferes 
constructively  with  the  previous  layer  reflections.  This  requires  that  reflections  from 
consecutive  layers  are  in  phase  with  respect  to  each  other  ( 2mn ;  m=0,l,2..  phase  difference 
between  reflected  waves).  Interference  between  reflected  waves  depends  upon  the  Optical 
Path  Difference  (A)  in  the  routes  taken  between  the  two  interfaces,  and  the  Fresnel  phase 
change  upon  reflection  from  an  interface.  Between  two  consecutive  layers,  A  is  given  by 
(from  Figure  8): 

A  =  n2\AB\  +  n2 \BC\  -  nx\AD\  ( nm )  (3) 

which  reduces  to 

A  =  ln2d  cos (Qt)  (nm)  (4) 

where  (Figure  8)  Qt  is  the  transmitted  angle  in  the  layer,  d  is  the  thickness  of  the  layer,  and 
nz  is  the  refractive  index  of  the  layer.  The  phase  difference,  512,  between  two  rays  from 

consecutive  layers  is  obtained  by  multiplying  A  by  the  wavenumber  k  and  then  adding  the 
Fresnel  phase  change  upon  reflection  from  each  layer. 


k\  +  <p]+(p2  =  — n2d  cos(0()  +  n  ( radians ) 

A 

(5) 

2  m7r  =  — n2d  cos(6?( )  +  n  ( radians ) 

A 

(6) 

4  TC 

2  (m  -  \)n  =  — n2d  cos(f9t )  ( radians ) 

(7) 

nx  sin  =  n2  sin  02  (unitless) 

Snell’s  Law  (8) 

2-2 


Figure  8.  Two  boundary  problem 


Table  1.  Fresnel  Phase  Change  Upon  Reflection 


S-Polarization  P-Polarization 


ni>nt 

0 

7t 

nt>ni 

Jt 

0 

Table  1  shows  the  Fresnel  phase  changes  applicable  to  the  RCLED  design.  A  more 
detailed  listing  is  available  in  OPTICS,  by  Hecht  [23:100],  The  Fresnel  phase  change, 
/S.(p-(p2  -<px,  sums  to  ±7t  (from  Table  1),  so  for  constructive  interference  kA  must  sum  to  a 

constant  phase  difference  of  7t.  Equation  6  shows  two  criteria  for  a  constant  phase  difference 
of  n  between  rays  reflected  from  different  layers.  First,  that  the  phase  difference  between 
two  consecutive  layers  ( 8n )  is  a  function  only  of  incident  angle  (through  Snell’s  law)  if  the 
layers  have  a  thickness  of  d  =  \/4n2.  Second,  that  the  design  wavelength  thickness  (A0/ 4) 

criteria  for  the  thin  layers  is  true  only  at  normal  incidence.  The  peak  reflected  wavelength 
(wavelength  required  for  2m /r phase  shift)  decreases  with  increasing  angle  0t.  Intuition 

might  suggest  that  as  the  pathlength  through  each  layer  increases  with  increasing  angle, 
the  design  wavelength  should  also  increase. 
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The  maximum  DBR  reflectance  occurs  at  the  design  wavelength.  A  simple  formula, 
neglecting  absorption,  to  determine  the  maximum  reflectance  for  normal  light  incidence 
(0  deg)  of  a  series  of  N  periodic  high-low  index  (HL)  X/2  layers  deposited  on  a  substrate  is 
given  by[24]: 


R»  = 


f  n  ^ 


2N\ 


nu 


\nL. 


1  + 


(  n  ^ 


f ..  \ 


IN 


n 


H 


\nLJ 


(unitless) 


(9) 


•  ns  is  the  real  refractive  index  of  the  substrate  layer 

•  no  is  the  real  refractive  index  of  the  incident  layer 

•  nH  is  the  real  refractive  index  of  the  high  index  XI 4  layer 

•  nt  is  the  real  refractive  index  of  the  low  index  XI 4  layer 

•  N  is  the  number  of  periods  of  HL  repeating  XI2  layers 

As  shown  for  an  example  Al^Ga^As  /Al^Ga^As  DBR  in  Figure  9,  the  reflectance 

increases  asymptotically  to  1  for  an  increasing  number  of  HL  layers.  The  reflectance  is  also 
larger,  for  a  given  number  of  layers,  for  larger  HL  index  difference,  i.e.  tsn-nH-nL.  Figure 

9  shows  that  high  reflectance  (R>0.99)  is  reached  with  about  25  HL  periodic  X/2  layers.  The 
stopband  bandwidth  of  the  reflectance  spectrum  may  be  calculated  from: 


AA  —  An 


-arcsin! 


n 


nH~nL 

\nH  +nL) 


(nm) 


(10) 


where  Xo  is  the  design  wavelength  of  the  DBR  [10].  For  the  example  of  Figure  10,  the 
bandwidth  should  be  83.15nm.  It  is  actually  ~95nm.  The  discrepancy  is  due  to  Equation  10 
giving  an  asymptotic  limit  for  an  increasing  number  of  HL  layers.  Figure  11  shows  that  at 


60  HL  layers  the  bandstop  has  reached  the  asymptotic  limit  of  Equation  10. 
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Figure  9.  Calculated  reflectance  vs.  number  of  HL  layers  for 
AlxHGaI  x[  As  /AlxLGal  xLAs  DBR  on  GaAs  substrate.  X0  =  850nm. 

An  =  nH-nL 


Figure  10.  Sample  calculated  reflectance  spectra  for  Al0 15Ga0  85As 
/AlAs  DBR  with  38  HL  periods  on  GaAs  substrate  (absorption 
effects  included)  XQ  =  850nm 
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Figure  11.  As  the  number  of  HL  half  wave  periods  of  Figure  10  are 
increased  the  bandstop  region  of  the  DBR  is  narrowed  to  the  limit 
of  Equation  10 


2.1.1  Compositional  Grading  of  DBRs.  In  AlGaAs  DBRs  most  of  the  resistance 
comes  from  the  hetero-interfaces  (between  the  HL  layers)  in  p-doped  mirrors  [13;  15:54-63]. 
Various  methods  of  decreasing  resistance  are  used.  One  method  of  reducing  resistance  is  to 
heavily  dope  the  p-mirror.  However  the  disadvantage  is  that  free  carrier  absorption 
increases,  and  this  decreases  quantum  efficiency.  The  method  used  in  this  thesis  is  to 
compositionally  grade  the  HL  interfaces  of  the  DBR.  For  a  compositionally-graded  DBR 
quarter  wave  layer  to  meet  the  Bragg  criteria  the  following  equation  must  be  true. 

|  «(z)Ao  dz  =  (nm)  (1 1) 

where  the  lower  limit,  zi,  is  the  start  of  that  particular  layer  of  growth  and  the  upper  limit, 
Z2,  is  the  end.  The  subscript  X0  denotes  that  this  relationship  is  true  only  at  the  design 
wavelength. 
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Compositionally  grading  the  DBR  decreases  resistance  because  abrupt  conduction  and 
valence  band  energy  offsets  at  hetero-junctions  act  as  barriers  to  current  flow  [15:54-55]. 
By  grading  the  interfaces  the  junctions  are  less  abrupt  and  resistance  to  current  flow  is 
decreased. 

Modern  deposition  systems  such  as  Metal-Organic  Vapor  Phase  Epitaxy  (MOVPE),  and 
Molecular  Beam  Epitaxy  (MBE)  create  semiconductor  structures  an  atomic  layer  at  a  time. 
It  is  therefore  possible  to  vary  the  composition  of  the  deposited  structure,  as  it  is  being 
built,  by  varying  the  relative  ratios  of  the  constituent  parts.  Figure  12  shows  example 
profdes  of  HL  layers  for  comparison  purposes.  Each  grading  type  has  the  same  value  of 
maximum  high  index  and  minimum  low  index.  The  only  difference  is  the  grading  scheme 
employed  between  these  index  upper  and  lower  limits. 


Figure  12.  Index  profile  for  periodically  graded  layers  of 

(A1XHGai-XH)0.5In0.5P/(A1xLGai.xL)0.5In0.5P  with  XH  =  01’  XL  =  0  9  and 

nH(A0)  =  3.429,  nL0y  =  3.171  where  Z0=  680nm 

Figure  13  shows  the  spectra  of  the  different  compositionally  graded  DBRs  of  Figure  12. 
The  material  system  is  (AlxGal  x)05In05P  with  xH  =  0.1  and  xl  =  0.9.  The  reflectance 

spectrum  of  a  DBR  depends  on  how  wide  the  high  and  low  index  composition  “plateaus”  in 
Figure  12  are  allowed  to  be.  Wider  “plateaus”  mean  a  steeper  grade,  and  higher  reflectance 
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values.  Narrower  “plateaus”  mean  a  less  steep  grade  and  lower  reflectance  values.  In 
Figure  13  it  is  easy  to  see  that  the  sinusoidal  and  step  grading  give  the  higher  maximum 
reflectances  because  they  have  wider  “plateaus”. 

Figure  14  shows  the  energy  bandgaps  and  lattice  constants  for  both  the  AlGaAs  and 
AlGalnP  families  of  III-V  compound  semiconductors.  The  lines  connecting  the  binary 
compounds  give  the  energy  bandgap  and  lattice  constant  for  the  ternary  compounds  made 
from  a  compositional  variation  of  the  elements  [15].  For  example,  for  the  system  A^Ga^As, 

at  x=  0,  it  has  the  properties  of  GaAs,  whereas  at  x  =  1,  it  has  the  properties  of  AlAs.  An 
important  feature  of  the  AlxGa1.xAs  compound  family  is  that  its  lattice  constant,  a,  hardly 

changes  as  jc  varies  from  0  to  1.  This  means  that  AlxGa1-xAs  DBRs  may  be  grown  with 
relative  impunity  with  regard  to  composition.  Lattice  matching  is  an  important  parameter 
in  growing  DBRs  because  growing  lattice  mismatched  layers  leads  to  lower  quality 
epitaxial  material  with  cracks,  defects,  and  dislocations.  If  growing  (AlxGa1.x)yInyP 
compounds  on  GaAs  substrates,  for  example,  more  care  must  be  taken  as  the  lattice 
constant  has  a  large  variation  with  x  and  y.  This  material  is  approximately  lattice  matched 
to  GaAs  when  y^0.5  and  0<x<l. 


2-8 


Energy  Bandgap  (eV) 


Wavelength  (ran) 


Figure  13.  Calculated  reflectance  spectra  from  50  graded 
(A^Ga^Q  5InP  HL  periods  with  a  GaAs  substrate  (n~3.7). 

xH  =  0.1,  xL  =  0.9  and  nH(/l0)  =  3.429,  nL(A0)  =  3.171  where 
X0  -  680nm 


Figure  14.  Energy  bandgap  vs.  lattice  constant. 


2.1.2  Complex  Refractive  Index.  Part  of  the  problem  of  matching  the  DBR  to  a 
design  wavelength  is  ensuring  that  the  DBR  does  not  attenuate  the  light  through 
absorption.  Absolution  occurs  when  the  photon  energy  is  greater  than  the  bandgap  energy 
of  the  DBR  layers.2 


When  the  photon  energy  is  greater  than  the  DBR  bandgap  energy,  photon  absorption 
must  be  considered.  The  complex  refractive  index  N  used  in  this  context  consists  of  the  real 
refractive  index  n  and  the  complex  extinction  coefficient  k  (N-n-iK).  The  form  of  the  electric 
field  wave  equation,  for  a  wave  traveling  in  the  x-direction  is  given  by: 


E  =  E0  exp<  /' 


cot - 


2n{n  -  itc) 


1  \ 


(Volts/m) 


(12) 


The  imaginary  part  of  the  complex  refractive  index,  k,  becomes  a  real  number  and 
denotes  an  attenuation  as  seen  in  Equation  13. 


£  =  £0exp(-^x}exp 


cot- 


2  n{ri) 


(Volts/m) 


(13) 


The  relationship  between  the  extinction  coefficient  k  and  the  Lambert-Beer  absorption 
coefficient  a  (m  ')  is  given  by  K-ca  !2co  where  c  represents  the  speed  of  light  in  a  vacuum, 
and  ©  is  the  angular  frequency  of  the  incoming  light. 


2  Photon  energy  and  wavelength  are  connected  by  E  —  he  /  X .  For  wavelength  given  in  fim  the 
numbers  simplify  to  E  =  1.24  /  X  (eV). 
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Figure  15  and  Figure  16  are  derived  from  work  by  Adachi  and  Aspnes.  The  refractive 
index  values  (real)  for  energies  less  than  1.5eV  come  from  extrapolation  equations  by 
Adachi.  The  complex  refractive  index  plot  is  based  on  data  presented  by  Aspnes.  Aspnes 
presented  his  data  numerically  for  varying  compositions  and  energy  values.  The 
compositions  varied  from  x  =  0.1  to  x  =  0.8.  The  energy  values  varied  from  1.5  eV  to  3.4  eV 
in  steps  of  0. 1  eV.  A  cubic  spline  function  was  used  to  interpolate  between  energy  values  for 
the  composite  Adachi/Aspnes  data  set .  A  linear  fit  was  used  to  interpolate  between  known 
compositional  values  when  the  composition  was  one  not  measured  by  Aspnes.  Figure  16 
shows  that  the  extinction  coefficient  vis  strongly  dependent  on  both  the  fractional 
composition  of  AlxGal  icAs  and  on  the  photon  energy. 


Figure  15.  Real  refractive  index  of  A^Ga^As 
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Figure  16.  Imaginary  refractive  index  of  AlxGalxAs 

When  absorptance  is  included,  the  form  of  the  phase  delay  changes  from  Equation  5. 

The  phase  delay  with  absorptance  is  given  by: 

Sn  =  — d2 (n\  -K22-nl  sin2  0O  -  2 in2K2)V2  ( radians )  (14) 

X 

where  n0  and  00  are  the  real  refractive  index  and  incidence  angle  of  the  incident  medium 

[10;20]. 

2.2  Fabry-Perot  Cavities 

A  Fabry-Perot  cavity  consists  of  two  parallel  planar  mirrors,  with  reflectances  R1  and 
R2.  The  classic  Fabry-Perot  cavity  transmits  wavelengths  that  are  defined  by  the 
separation  between  the  plates  (dcav).  The  condition  is  dcm  =  mX0  / (2 ncav)  for  transmission 
through  the  cavity,  where  m  is  the  longitudinal  mode  number,  Xo  is  the  design  wavelength, 
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and  ncau  is  the  refractive  index  inside  the  cavity  [25].  The  transmissivity  for  these  resonant 
wavelengths  is  very  high  (~1)  if  Rs  ~R2~  1.  The  resonant  condition  is  shown  by  Figure  17 

where  increasing  reflectances  on  the  mirrors  create  a  higher  Q  (quality  factor)  cavity.  The 
formula  that  describes  the  transmittance  of  a  Fabry-Perot  cavity  with  R}  =  R2  =  R  is: 


T  = - v — — — o -  (unitless) 

(1  -  R)2  +  4R  sin2  (<p) 


where  tp  as  a  function  of  the  angle,  0cav,  inside  the  cavity  is  given  by 


<P  =  —ncavdcav  cos  8^  (radians) 


<j>(in  multiples  of  ti) 


Figure  17.  Transmittance  characteristics  of  Fabry-Perot  cavity 
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The  free  spectral  range  (FSR)  is  the  distance  in  nm  between  the  longitudinally 
propagating  modes.  The  FSR  is  only  affected  by  the  mirror  separation,  ncavdcav.  The  full- 

width  at  half-maximum  (FWHM)  is  defined  as  the  wavelength  spread,  AX,  about  the  design 
wavelength.  The  FWHM  decreases  for  increasing  reflectance  of  R;  and  R2.  The  cavity  Q  is 
the  ratio  of  the  FWHM  to  the  design  wavelength  X ^  Higher  Q  values  mean  greater  mirror 
reflectances  and  narrower  Fabry-Perot  cavity  transmission  passbands  as  seen  in  Figure  17. 

2.2.1  RCLED  Microcavity.  The  RCLED  Fabry-Perot  cavity  is  formed  by  growing  a 
DBR  on  a  GaAs  substrate,  followed  by  an  active  optical  cavity  of  an  integral  number  of  half 
Bragg  wavelengths,  and  then  another  DBR.  This  is  illustrated  in  Figure  18  where  metal  • 
contacts  with  emission  apertures  have  been  placed  on  the  upper  DBR  and  a  blanket  metal 
has  been  put  on  the  substrate  to  form  an  excitation  current  path.  The  cavity  of  the  RCLED 
is  referred  to  as  a  microcavity  because  its  small  size  (~Ao)  changes  some  of  the  characteristic 
Fabry-Perot  parameters.  [7;  8;  26;  27], 

The  resonant  condition  for  a  microcavity  with  DBRs  is  modified  by  the  penetration 
depth  of  the  electric  field  intensity  into  the  DBR  mirrors.  The  penetration  depth  is  a 
parameter  used  to  replace  a  distributed  mirror  structure  having  linear  phase  variation  for 
changing  frequency  with  a  fixed  phase  mirror  (Figure  19). 
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Fabricated  Devices 


Ti/Au  pcontact 


Upper 
DBR  1 

Cavit^ 

i 

Lower 

DBR 


Ge/Au/Ni/Au  n  contact 

Figure  18.  RCLED  structure  showing  microcavity,  DBRs,  and 
devices  (not  to  scale) 


Figure  19.  Penetration  depth  in  distributed  mirrors 
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The  condition  for  the  field  reflected  from  the  lower  distributed  mirror  and  the  lower 
fixed  mirror  to  have  equal  phase  is  given  by  [26]: 


(a-(O0)Tr 


(radians) 


(17) 


given: 


0O  =  co0rr  (radians) 

(IB) 

=  ne„tanh(KQ;  (radiam) 

(19) 

C  K 

3  =  2  nc  /  XQ  (radians /sec) 

(26) 

•  k  :  coupling  coefficient 

•  r  :  DBR  reflection  delay 

•  n  :  effective  refractive  index  of  cavity 

cav 

•  l :  thickness  of  DBR  structure 

•  0O  :  the  phase  change  at  the  fixed-phase  mirror. 

The  resonant  condition  for  the  Fabry-Perot  microcavity  is  that  the  wavefront  phase  be 
self  consistent  (equal)  for  a  round  trip  in  the  cavity: 


Y  >  “  0'—  "  =  2mn  (radians)  (21) 

dcff  =  ncmdcm  +nupperdupper  +nlowerdlower  ( nm )  (22) 

where  dupper  and  dlower  represent  the  electric  field  intensity  penetration  depths  into  the 

upper  and  lower  DBRs  respectively.  Equation  21  suggests  that  the  new  effective  length  of 
the  microcavity  should  change  the  primary  resonance  condition  (design  wavelength). 
However,  for  a  microcavity  that  is  wavelength  matched  with  the  DBR  (as  in  an  RCLED) 
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this  is  not  the  case.  The  DBR  always  reflects  it’s  design  wavelength  with  either  0  or  n  phase 
change  and  this  causes  Equation  21  to  give  dcav  =  deff  at  the  design  wavelength  at  normal 

incidence. 

The  characteristic  narrowness  of  the  spectral  output  of  the  RCLED  is  caused  by  the 
microcavity.  Only  certain  modes  propagate,  and  only  those  modes  may  be  excited  by  the 
gain  spectrum  of  the  quantum  wells.  Therefore,  it  is  essential  to  know  the  mode  structure 
of  the  RCLED  microcavity  [9]. 

2.2.2  RCLED  Microcavity  Parameters.  The  mode  structure  of  the  microcavity  is 
defined  by  it’s  characteristic  parameters.  That  is,  the  quality  Q,  the  Free  Spectral  Range  ■ 
FSR,  the  Finesse  F,  and  the  photon  lifetime  rp  in  the  cavity.  Because  of  the  distributed 

nature  of  DBRs,  these  characteristic  parameters  change  for  RCLEDs. 

Q  =  —  ( unitless )  (23) 

AA 

FSR  =  — — —  =  A2°  (nm)  (24) 

2nd  _  2 nd,ff 


F  = 


FSR 


FWHM  2ncavdeffAA 


(unitless) 


(25) 


QK  _  2ncavdejSr 


2nc  cil-R^Rz) 


(sec) 


(26) 


2.3  Reflectance  Calculations  for  RCLED 

The  Mathcad  5.0+®  document  used  to  calculate  the  reflectance  spectra  for  this  thesis 
uses  a  finite  element  approach  to  calculating  device  reflectance.  Each  quarter-wave  (A/4) 
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layer  is  broken  into  N#uW  separate  sublayers.  More  sublayers  give  greater  accuracy  in 
modeling  the  smooth  compositional  grading  of  the  X/4  layers.  However,  increasing  NsuWayer 

results  in  increased  processing  time.  For  this  reason,  the  program  breaks  graded  areas  into 
sublayers,  but  treats  the  “plateaus”  in  the  flayers  as  a  single  sublayer  (see  Figure  20). 
The  reflectance  program  using  the  Macleod  matrix  method  makes  the  following 
assumptions  about  the  multilayer  device  being  modeled: 

1.  Non-magnetic  materials 

2.  Smooth  interfaces  without  scatter 

3.  Plane  parallel  layers  with  infinite  lateral  extent 

4.  Semi-infinite  incident  (air)  and  substrate  (GaAs)  media 

5.  Linear  wave  equation. 

The  program  can  handle  arbitrary  compositional  grading  with  complex  refractive  index 
information.  The  condition  for  setting  the  thickness  of  the  Nsubiayer  separate  sublayers  is  a 
discrete  form  of  Equation  11. 

N sublayer  1 

2X(z)^  =  -r  ( nm )  <27) 

r=  1  ^ 

where  nr(z)  denotes  the  compositional  grading  in  the  growth  direction,  dr  is  the  sublayer 
physical  thickness,  and  the  subscript  X0  means  the  equation  should  be  calculated  at  the 
design  wavelength,  X0.  This  ensures  that  each  A/4  layer  is  the  correct  optical  thickness  at 
the  design  wavelength.  The  criteria  for  the  optical  cavity  of  the  RCLED  is 

^ cavity 

=  (nm)  (28) 

r=l  ^ 


2-18 


where  Ncavity  is  the  number  of  sublayers  modeled  for  the  cavity  and  NUM  is  the  width  of  the 
cavity  in  half- wave  lengths. 


The  reflectance  program  listing  is  included  in  Appendix  B. 
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Figure  20.  Scheme  for  dividing  HL  layers  into  sublayers  in 
reflectance  program 


2.3.1  Example  Calculation.  Figure  21  shows  the  calculated  reflectance  pattern  for  an 
RCLED  with  38  periods  in  the  lower  DBR,  a  cavity  of  1  Bragg  wavelength,  and  6  periods  in 
the  upper  DBR  (Figure  18).  The  HL  periods  are  composed  of  Al0 15Ga0  85As/ALAs  with 
interfaces  that  are  linearly  graded  over  12nm.  The  reflectance  of  the  upper  DBR  is  less 
than  that  of  the  lower  DBR  in  order  to  favor  emission  through  the  upper  surface.  The  GaAs 
substrate  is  absorbing  for  wavelengths  shorter  than 

X  <  ~  870nm.  As  expected,  the  characteristic  Fabry-Perot  dip  is  seen  to  move  toward 
decreasing  wavelengths  with  increasing  incidence  angle. 
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Reflectance 


Figure  21.  Reflectance  spectrum  of  RCLED  with  6  HL  periods  of 
Al0 15Ga085As/AlAs  in  the  upper  DBR,  38  HL  periods  of 

Al0 15Ga0  85As/A1As  in  the  lower  DBR  and  a  1A,,  at  A,0  =  850nm, 
cavity  (Calculated  for  TE  polarization) 

Compare  the  characteristics  of  a  RCLED  with  a  I  a  cavity  at  850nm  with  ricav  =  3  and 
Rapper  =  0.5,  Riower  =  0.99  to  a  typical  He-Ne  laser  cavity  with  a  20cm  cavity  at  X  -  632.8nm, 
ricav  =  1,  Rapper  =  0.95,  and  Riower  =  0.99.  The  results  are  given  in  Table  2. 


Table  2.  Comparison  of  Selected  Cavity  Parameters. 


FSR 

Finesse 

Photon  Lifetime 

RCLED  Cavity 

425  nm 

8.889 

11.22  femto  sec 

HeNe  Cavity 

0.001  nm 

102.422 

22.41  nano  sec 
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A  major  advantage  of  the  RCLED  is  the  width  between  the  modes  (FSR).  For  the  He-Ne 
laser  at  632.80025nm  the  mode  number,  q,  is  632,111  and  at  632.80125nm  q  is  632,110. 
Therefore  the  gain  spectrum  is  wide  enough  for  many  modes  to  propagate  and  multimode 
operation  is  the  norm.  For  the  RCLED  at  850nm  the  mode  number  q  is  2.  For  q  -S  the 
wavelength  is  566nm.  Thus  the  RCLED  is  single  mode  because  the  microcavity  formed  by 
DBRs  is  wavelength  dependent  and  does  not  reflect  sufficiently  at  566nm,  and  the 
quantum  well  gain  spectrum  does  not  extend  to  566nm. 

2.4  Quantum  Well  Layers 

The  last  wavelength  selection  capability  occurs  inside  the  quantum  well  (QW)  layers.  • 
Saleh  and  Teich  describe  a  layer  as  being  quantum  in  nature  when  it  is  “smaller  than  the 
de  Broglie  wavelength  of  a  thermalized  electron  [28].”  For  the  present  study,  quantum  well 
layers  must  be  25nm  thick,  or  less. 

The  wavelength  of  the  light  emitted  from  a  semiconductor  diode  is  dependent  upon  its 
energy  bandgap.  Quantum  wells  have  discrete  allowed  carrier  (electron  and  hole)  energy 
levels  in  the  direction  of  growth.  These  discrete  levels  mandate  that  electron-hole 
recombinations  are  only  allowed  for  certain  energies.  Therefore  the  art  of  bandgap 
engineering  (quantum  tailoring)  involves  the  exploitation  of  individual  layers  of  atoms  to 
realize  properties  that  would  not  exist  otherwise  [29;30]. 

2.4.1  Ternary  and  Quaternary  Compound  Materials.  Ternary  (three  element)  and 
quaternary  (four  element)  compound  materials  are  used  by  device  makers  to  realize  desired 
energy  bandgaps  and  lattice  constants  [15;  29]. 

Figure  22  shows  the  range  of  useful  wavelengths  that  can  be  achieved  by  some 
compound  materials.  The  object  of  interest  to  this  thesis  is  the  AlxGal  xAs  family,  capable  of 
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emitting  light  from  0. 7 /urn  to  0.9 /^m  wavelength  depending  on  the  value  of  x  (0  <  x  <  0.45). 
The  composition  x  is  cutoff  at  0.45  because  AlxGalxAs  becomes  an  indirect  semiconductor 

What  makes  these  composite  compounds  so  useful  is  that  the  energy  bandgap  may  be 
changed  in  situ  by  varying  the  composition  with  time  during  growth. 

2.4.2  Allowed  Energy  Levels.  Figure  23  schematically  shows  the  approximately 
parabolic  shells  that  define  the  allowed  energies  of  the  quantum  well  electrons  and  holes 
close  to  the  subband  energy  levels  [28]. 


Wavelength  (jam) 


Figure  22.  Range  of  possible  emission  wavelengths  for  selected 
compound  semiconductors 
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(a)  (b)  *  (c) 


Figure  23.  (a)  Geometry  of  quantum  well,  (b)  Energy  level 
diagram  for  electrons  and  holes  in  quantum  well,  (c)  Cross  section 
of  E-k  relation 

In  bulk  material  the  E-k  (energy-momentum)  relationships  for  electrons  near  the 
bottom  of  the  conduction  band  and  holes  near  the  top  of  the  valence  band  are  given  by 


r  r  n2k2 

E-E.-v - 

(eV) 

(29) 

2mc 

_  IT  ***' 

E-Ev - 

(eV) 

(30) 

2  mv 

k2  =k2  +k2  +k2 

(m-2) 

(31) 

where  m  and  m  are  the  effective  masses  in  the  conduction  band  and  valence  band 

c  v 

respectively,  and  the  subscripts  {x,y,z}  denote  a  coordinate  system  choice.  In  a  quantum 
well  the  energy  levels  become  discrete  in  the  direction  of  growth,  but  remain  continuous  for 
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the  transverse  (x,y)  directions  (perpendicular  to  growth).  The  E-k  relationships  for  the 
conduction  and  valence  bands  for  the  quantum  well  grown  in  the  z  direction  are: 


E1=E.+Ek+^ i-  (eV)  (32) 

2  mc 

%2k2 

=EV-E  (eV)  (33) 

2  mv 

K=K+K  C m -*)  (34) 

Therefore  all  electron-hole  recombinations  leading  to  photon  emission  are  at  energies 
defined  by: 

E^=E2-E,=E„v+E,u+E,„  +  ?£E  (eV)  (35) 

where  Egap  is  the  natural  bandgap  of  the  quantum  well  material  and  the  reduced  mass  rnr 
is  calculated  using  l/mr  =  \/mc  +\/mv.  From  Figure  23,  electron-hole  recombinations, 

which  lead  to  photon  emission,  are  not  allowed  at  the  bandgap  energy  in  a  quantum  well, 
but  must  take  place  at  a  higher  energy.  The  first  allowed  energy  transition  is  between  qcz  = 

1  in  the  conduction  band  and  qvz  =  1  in  the  valence  band  (q(c-v)z  =  1).  With  careful  selection 
of  materials  and  QW  geometry  it  is  possible  to  choose  emitted  wavelength.  Figure  24  shows 
the  wavelength  dependence  versus  well  width  of  the  q(c-v)z  =  1  (electron  -  heavy  hole) 
transition  for  a  AlxGal  xAs  quantum  well  with  a  Al0 15Ga085As  barrier  region. 
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Figure  24.  Dependence  of  emission  wavelength  on  AlxGalxAs  well 
width  and  aluminum  content  x  with  an  Al0 15Ga0  85As  barrier 
region 

Figure  25  and  Figure  26  show  the  wavefunctions  describing  the  location  of  the  carrier  in 
the  conduction  and  valence  bands  of  the  quantum  well.  The  figures  are  calculated  by 
solving  the  one -dimensional  time-independent  Schrodinger  equation  [19]  for  a  15nm  GaAs 
quantum  well  with  Al03Ga0  7As  barrier  layers  (Appendix  C).  There  are  2  allowed  electron 

levels,  7  allowed  heavy  hole  levels  (4  are  shown)  and  3  allowed  light  hole  levels  (none  are 
shown).  The  probability  of  a  transition  between  levels  is  proportional  to  the  correlation  of 
the  energy  level  wave  functions  given  by 

Ptrans  00  J ^ m(Z)'V ( UnitleSS )  (36) 
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Figure  25.  Wavefunctions  of  energy  levels  in  conduction  band 
calculated  for  a  15nm  GaAs  quantum  well  with  Al03Ga07As 

barrier  layers 
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Figure  26.  Wavefunctions  of  heavy  hole  energy  levels  in  valence 
band  calculated  for  a  15nm  GaAs  quantum  well  with  Al03Ga0  7As 

barrier  layers 
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The  wave  functions  are  given  by  'F  oc  cosAz  for  odd  energy  levels  (1,  2,  3...)  and 
Y  oc  sinAz  for  even  energy  levels  (2,4,6...)  where  A  is  the  wave  number  of  the  energy  level 
and  z  is  the  distance  across  the  quantum  well.  Thus  transitions  between  levels  m  and  n  are 
more  likely  for  m  =  n  than  transitions  between  levels  where  m^n .  For  the  purpose  of  this 
thesis,  rigorous  transition  rules  are  applied.  Ptrans  is  assumed  to  be  zero  for  m^n  and  only 
electron-heavy  hole  recombinations  are  considered  .  Electron-light  hole  recombinations  are 
ignored  in  this  thesis  because  they  occur  at  different,  higher  energies  and  therefore  will  not 
be  enhanced  by  the  microcavity. 

2.4.3  Quantum  Well  Gain  Spectrum.  The  optical  gain  available  to  a  photon  of  energy 
E  in  a  quantum  well  is  given  by  [28]: 

g(E)  =  /-p(E)fs(.E)  <m-0  (37) 

8^rrr 

•  p(E)  is  the  2D  optical  joint  density  of  states 

•  E  is  the  photon  energy 

•  fg(E )  is  the  fermi  inversion  factor 

•  Tr  is  the  electron-hole  recombination  lifetime. 


For  a  quantum  well  with  i  electron  energy  levels  and;  hole  energy  levels  Equation  37 
becomes[16;31]. 


g(E ) 


q2\Mf 


Es0mlc0hnL. 


<cmi) 


(38) 


» ,j 


where 

•  i,  j  are  integers  referring  to  the  quantized  electron,  heavy  hole,  and  light  hole  states 
respectively 

•  E  is  the  photon  energy 
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.  |  Mp  is  the  transition  matrix  element  [31] 

•  m  is  the  reduced  mass  of  the  i,j  transition 

r  1>J 

•  C  is  the  strength,  or  degree  of  wavefunction  correlation,  of  the  i,j  transition 

hj 

•  A  .  accounts  for  the  polarization  dependence  of  the  ij  transition 

V 

•  f  and  /v  are  the  Fermi-Dirac  distribution  functions  for  the  conduction  and  valence 
bands 

•  <P  is  the  Heaviside  step  function  {<b(x)  =  0,x  <  1  and  l,x  >  0}. 

•  n  is  the  effective  group  refractive  index  of  the  quantum  well 

•  Lz  is  the  thickness  of  the  quantum  well. 

In  order  to  include  the  spectral  broadening  of  each  recombination  transition  due  to 
interactions  with  phonons  and  other  electrons,  the  gain  expression  (Equation  38)  is 
convolved  with  the  Lorentzian  spectral  lineshape  function 


Lift* 


1 

7t(E-h(0)2  +(h/Tin)2 


(eV-0 


(39) 


to  get  the  final  calculated  gain 

oo 

G(hco)  =  |  g(ho)  -  Etrans)L{EtransyEtrans  (cm1)  (40) 

-oo 

The  optical  gain  for  the  RCLED  depends  on  the  quasi  fermi  levels,  Efc  and  Efv,  which  are 

a  function  of  carrier  density  in  the  well.  The  relationship  between  the  carrier  densities; 

N electrons,  Phoies,  and  the  energy  difference  between  the  ith  quantized  electron/hole  energy  level 

and  the  quasi- fermi  levels,  assuming  parabolic  E-k  bands,  is  given  by: 

^ electrons  =  ^{1+eXpH^  ~Efc)lkJ )}  (cm*)  (41) 
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(42) 


Ajrk  T 

Pholes  =  -777—  Z  ln(1  +  eXP(-(^v,  -Efi,)!  kBT)}  (cm-5) 

where  4  is  the  width  of  the  quantum  well.  For  this  thesis  charge  neutrality  is  assumed  and 

N electrons  =  Phoies.  Given  N, electrons  and  Phoies,  Equations  41  and  42  may  be  solved  for  the  quasi  - 
fermi  levels.  E,  and  Efy  are  then  put  back  into  the  gain  equation  and  used  to  determine 

G(f)co). 

Figure  27  shows  the  calculated  gain  for  an  8nm  GaAs  quantum  well  with  Al0 15Ga085As 
barrier  layers  (same  parameters  as  Figure  25).  As  stated  before  only  the  q(c-v)z=l  electron 
to  heavy  hole  transition  is  being  considered  as  contributing  to  the  gain. 

The  width  of  the  gain  spectrum  is  defined  as  the  energy  difference  between  the  q(c-v)z=l 

electron  to  heavy  hole  transition  and  the  separation  of  the  quasi-fermi  levels.  Put  another 
way,  the  quantum  wells  provide  gain  for  photon  energies  in  the  range 

Ex<E<Efc-Ep  (elO  (43) 

where  E}  is  the  transition  energy  for  q(c-v)z  =  1. 

2.4.4  Quantum  Well  Rate  of  Spontaneous  Emission.  The  spectral  rate  of 
spontaneous  emission  from  the  quantum  well  is  the  summation  over  all  possible  energies  of 
the  product  of  the  transition  probability  density,  the  emission  probability,  and  the  2D 
optical  density  of  states  that  can  interact  with  the  photon  [28]. 

R(E)  =  —  p(E)fe(E)  {J1  cm3 sec1)  (44) 

G 

•  p(E)  is  the  2D  optical  joint  density  of  states 

•  E  is  the  photon  energy 
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•  fe(E)  is  the  probability  of  photon  emission 

•  tr  is  the  electron-hole  recombination  lifetime. 

The  rate  of  spontaneous  emission  from  a  quantum  well,  in  accordance  with  Equation 
44,  is  given  by 


R{E)  = 


2nE\M\ ’■ 


£om0c0h  L, 


-llmnJCi.jAijfXl-f^{E-Euj)  (J^cm^sec1)  (45) 


Figure  27.  Gain  for  8nm  GaAs  quantum  well  with  a  Al0 15Ga085As 

barrier  layers  for  various  well  carrier  concentrations  showing  both 
calculated  (sharp)  and  convolved  (smoothed)  spectra  at  T=300K 


To  obtain  the  spontaneous  emission  at  a  given  photon  energy  Equation  45  must  be 
convolved  with  the  Lorentzian  line-shape  function  as  before.  The  quasi  fermi  levels  are 
calculated  as  before.  In  the  RCLED  microcavity,  it  is  the  spontaneous  rate  of  emission  that 
is  modified  by  the  microcavity.  This  can  also  be  thought  of  as  modification  of  the 
spontaneous  emission  lifetime  by  the  microcavity. 
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Figure  28.  Spontaneous  emission  rates  for  a  8nm  GaAs  quantum 
well  with  Al0 15Ga085As  barrier  layers  for  various  well  carrier 

concentrations  showing  both  calculated  (sharp)  and  convolved 
(smoothed)  spectra  at  T=300K 


2.5  The  RCLED  Structure 

The  Fabry-Perot  microcavity  and  the  quantum  well  layers  work  together  to  select  the 
emitted  wavelength  of  the  RCLED.  The  next  step  is  to  examine  some  parameters  that 
affect  the  RCLED. 

2.5.1  Spontaneous  Emission  Enhancement  Factor  E.  The  spontaneous  emission 
enhancement  factor  ,E,  of  a  microcavity  is  the  ratio  of  the  spectral  output  of  a  dipole 
emitting  into  a  microcavity  to  a  dipole  emitting  into  a  bulk  medium.  E  describes  the 
modification  of  the  emission  due  to  the  longitudinal  mode  structure  of  the  microcavity.  An 
electron-hole  recombination  may  be  considered  as  the  source  of  coherent  wave-packets 
emitted  in  opposite  directions.  These  wave  packets  evolve  independently  upon  reflection 
leading  to  interference  outside  the  microcavity.  The  complex  reflectivity  of  the  electric  field 
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at  a  DBR  is  quantified  by  the  amplitude  modifier  p  and  the  phase  change  0,  where  0<p<l 
and  O<0<27t.  The  incident  and  reflected  E-fields  are  related  by: 


E 


reflected 


=  Pei<>  Cadent 


(46) 


E  is  derived  in  Chapter  3  but  is  shown  here  for  completeness.  E  as  a  function  of 
wavelength  at  0°  incidence  is  given  by: 


E 


emitted 


(1  -  R,  )f  1  +  R,  +  ncav  (d  -  dx )  +  0 


qw\ 


1 +RA~  cos(  +  01+0. 


(47) 


•  Rl  is  the  reflectance  of  the  upper  DBR 

•  R2  is  the  reflectance  of  the  lower  DBR 

•  n  is  the  refractive  index  of  the  cavity 

cav 

•  d  is  the  width  of  the  cavity 

•  c?1  is  the  distance  between  the  quantum  wells  and  the  upper  DBR 

•  0X  is  the  phase  change  upon  reflection  by  Ri 

•  02  is  the  phase  change  upon  reflection  from  R2 

•  |  Eqw  | 2  is  the  intensity  of  the  E-field  emitted  by  the  quantum  well  layers. 

•  |  Eemitted  1 2  is  the  intensity  of  the  E-field  emitted  from  the  RCLED. 


d 


Air 


pit  0}  Upper  DBR 


d, 

Microcavity 

d2 


Lower  DBR 
GaAs  Substrate 


Quantum  Well  Layer  = 

Figure  29.  Graphic  representation  of  architecture  for  E  in 
Equation  47,  only  normal  propagation  is  considered  here 
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2.5.2  Electric  Field  Intensity  in  the  Cavity.  The  electric  field  confinement  factor  is  a 
measure  of  how  much  the  electric  field  overlaps  the  quantum  well  layers.  The  figure  of 
merit  for  vertical  structures,  such  as  an  RCLED,  is  the  longitudinal  confinement  factor 
(along  the  direction  of  growth)  Tz.  Fz  is  given  by: 

| E2(z)dz 

p  -  unitless )  (48) 

■=\E\z)dz  ( 

d 

where  d  is  equal  to  the  cavity  width  and  Active  denotes  the  extent  of  the  quantum  wells. 
Thus  the  comparison  is  made  with  the  effective  coupling  between  the  energy  in  the  cavity 
and  the  energy  in  the  quantum  wells.  Higher  confinement  factors  mean  greater  efficiency' 
and  lower  losses  in  the  cavity.  Figure  30  and  Figure  31  show  the  calculated  electric  field 
intensity  in  a  RCLED  with  a  IX  optical  cavity.  The  RCLED  is  designed  for  operation  at 
~850nm.  Clearly  the  RCLED  will  not  work  well  at  825nm. 


4  4.2  4.4  4.6  4.8  5  5.2  5.4  5.6  5.8  6 

Distance  from  substrate  (|4m) 


Figure  30.  Normalized  electric  field  intensity  inside  RCLED 
plotted  against  distance  from  substrate  for  X  =  841nm 
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Distance  from  substrate  (|um) 

Figure  31.  Normalized  electric  field  intensity  inside  RCLED 
plotted  for  against  distance  from  substrate  X  -  825  nm 

2.5.3  RCLED  Losses.  For  stimulated  emission  to  occur,  gain  in  the  cavity  must 
overcome  losses  due  to  absorption,  mirror  transmittance  and  other  sources.  The  threshold 
gain  necessary  for  lasing  in  vertical  devices  is  given  by  [32]: 

gJAca,  =atL  +  amLm  +  aoc(doc  -L)-  An(RA)  ( unitless )  (49) 

a,  =  aaYz  +  ascat  +(1-T ,)ap  (cm1)  (50) 

•  Fz  is  the  longitudinal  confinement  factor 

•  gth  is  the  threshold  gain 

•  dcav  is  the  width  of  the  cavity 

•  L  is  the  gain  length  (sum  of  width  of  quantum  wells) 

•  am  is  the  scattering/absorption  loss  within  the  DBRs 

•  Lm  is  the  penetration  depth  into  the  DBRs 

•  acc(doc-L)  accounts  for  absorptive  losses  in  the  passive  region  of  the  cavity 

•  aa  is  the  free  carrier  loss  within  the  quantum  wells 
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•  aP  is  the  free  carrier  loss  within  the  passive  layers 

•  ascat  is  scattering  loss  throughout  the  cavity  due  to  sidewall  roughness 

•  ln(RiRz)  accounts  for  losses  due  to  mirror  transmittance 

Given  that  the  only  major  difference  between  the  RCLED  and  the  VCSEL  is  the  upper 
mirror  reflectance,  IniRJUJ  is  the  most  important  term  in  comparing  these  similar  devices. 
An  example  is  an  RCLED  and  a  VCSEL  each  with  lower  mirror  reflectance  R2  ~  0.999.  The 
VCSEL  has  upper  mirror  reflectance  =  0.95  while  the  RCLED  has  Rj  =  0.5.  Both  have 
the  same  microcavity  architecture  and  the  same  number  of  quantum  wells.  Therefore  in 
Equation  49  all  cavity  and  scattering  losses  are  presumed  equal,  as  is  the  longitudinal 
confinement  factor.  This  reduces  Equation  49  to: 

— rcled)  ( unitless )  (51) 

Sth  RCLED 

Using  the  numbers  in  the  scenario  the  RCLED  has  3.116  times  the  threshold  gain 
necessary  for  lasing  as  the  VCSEL.  This  implies  that  the  RCLED  will  lase  only  if  the  upper 
DBR  reflectance  is  increased. 

2.6  Summary 

Chapter  2  presents  all  of  the  basic  theory  behind  the  Resonant  Cavity  Light-Emitting 
Diode  in  a  piecewise  fashion.  The  physics  of  the  basic  RCLED  components  are  examined 
separately.  In  a  real  device  all  of  these  components  must  be  designed,  and  built,  to  work 
together  at  the  design  wavelength.  Component  integration  is  the  focus  of  Chapters  3  and  5. 
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3.  Designing  The  RCLED 


3.1  Introduction 

In  order  to  examine  and  quantify  the  effects  of  a  microcavity  on  the  spectral 
emission  properties  of  a  light-emitting  diode  the  individual  components  of  the 
RCLED  must  be  examined  separately.  Two  devices  were  designed  for  this  purpose. 
The  first  device  consisted  of  a  microcavity  deposited  on  a  GaAs  substrate.  This 
device  gave  the  unenhanced  emission  spectrum  of  a  microcavity  with  four  quantum 
wells,  Eqw.  The  second  device  has  complete  upper  and  lower  DBRs  and  is  a  true 

RCLED.  The  emission  spectrum,  Eemitted,  is  enhanced  by  the  microcavity  and  the 
location  of  the  quantum  wells. 

By  comparing  emission  spectra  from  each  device  the  effects  of  enhancement  of 
the  microcavity  formed  by  the  two  DBRs  may  be  examined.  The  goal  of  this  chapter 
is  to  quantify  that  change  in  a  meaningful  manner.  Both  devices  are  shown  in 
Figure  32. 

3.2  Structure  1:  The  Cavity 

Figure  33  shows  the  as-grown  composition  profile  for  Structure  1.  The  cavity 
contains  four  8nm  wide  GaAs  quantum  wells  separated  by  lOnm  Al015Ga0  85As 

barrier  layers.  The  barrier  layers  are  wide  enough  that  negligible  coupling  occurs 
between  quantum  wells.  The  absorptive  and  gain  effects  of  the  quantum  wells  are 
assumed  to  be  four  times  that  of  a  single  one  (i.e.  uniform  carrier  injection  in  every 
quantum  well). 
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Note  that  Structure  1  has  a  single  quarter-wave  layer  above  and  below  the 
cavity.  This  is  in  order  to  properly  define  a  separate  confinement  heterostructure  for 
the  quantum  well  active  region.  The  energy  levels  of  a  single  8nm  GaAs  quantum 
well  surrounded  by  Al0 15Ga085As  barrier  layers  are  shown  in  Figure  34.  For  the 

purposes  of  this  thesis  only  interactions  between  the  electron  level  and  the  first 
allowed  heavy  hole  level  are  considered  (i.e.  e\  — »  hhl).  Figure  34  was  generated 
using  the  material  parameters  in  Table  3,  where  m0  is  the  mass  of  an  electron.  Etrans, 

the  calculated  transition  energy  for  Figure  34,  is  1.458  eV  or  850nm.  Figure  35 
shows  the  calculated  gain  spectrum  for  the  four  quantum  wells  as  a  function  of 
wavelength  for  various  carrier  densities.  Figure  36  shows  the  calculated  rate  of 
spontaneous  emission  versus  energy  for  various  carrier  densities,  assuming  uniform 
charge  distribution  in  the  quantum  well.  Figure  36  represents  the  expected 
unenhanced  spontaneous  emission  from  the  microcavity.  The  smoother  curve  in 
Figure  35  and  Figure  36  has  been  convolved  with  the  Lorentzian  function  to 
approximate  intra-band  scattering  effects. 


dcav  ^ 


3>.  Buffer 
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Structure  1 
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Figure  32.  Schematic  of  RCLED  structures 
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Figure  33.  Composition  diagram  of  Structure  1 


Etrans~  Allowed  Transition 


E  f  Energy  Bandgap  For  GaAs 

Eg2=  Energy  Bandgap 
For  Al0il5Gao_85As 


hh  =  Heavy  Hole  Energy  Level 

lh  =  Light  Hole  Energy  Level 
e  =  Electron  energy  level 
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Figure  34.  Calculated  allowed  energy  levels  for  quantum 
wells  of  Structure  1  at  T=300K 
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Table  3.  Material  parameters  for  Figure  34 


in  well 

in  barrier 

Mass  of  electron 

0.067mo 

0.079m0 

Mass  of  heavy  hole 

0.62m0 

0.641m0 

Mass  of  light  hole 

0.087m0 

0.096mo 

Energy  Levels 

eV 

AEC 

0.122 

AEv 

0.075 

Egl 

1.424 

Eg2 

1.536 

7  x  1018cm'3 
carriers  in  the  well 


5  x  1018cm'3 


3  x  1018cm"3 


750  760  770  780  790  800  810  820  830  840  850  860  870  880  890  900  910  920  930  940  950 

Wavelength  (nm) 

Figure  35.  Calculated  gain  of  an  8nm  GaAs  quantum  well 
surrounded  by  lOnm  Al0 15Ga0  8BAs  barrier  layers 


Wavelength  (nm) 


Figure  36.  Calculated  rate  of  spontaneous  emission  of  an 
8nm  GaAs  quantum  well  surrounded  by  lOnm 
AIq  15Ga0  85As  barrier  layers 

Figure  35  and  Figure  36  were  calculated  with  Mathcad  using  the  document  in 
Appendix  C  and  the  equations  presented  in  Chapter  2.4. 

3.3  Structure  2:  The  Cavity  with  Upper  and  Lower  DBRs  (RCLED) 

In  Structure  2  the  resonant  cavity  is  fully  formed  with  the  inclusion  of  upper  and 
lower  linearly  graded  DBRs.  The  composition  diagram  of  Structure  2  is  shown  in 
Figure  37.  Figure  38  shows  the  calculated  normal  incidence  reflectance  and  phase 
(with  absorption)  for  the  RCLED  structure  in  Figure  37.  An  interesting  feature  is 
the  blue-shifting  of  the  Fabry-Perot  dip  predicted  by  the  model  from  the  as-designed 
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wavelength  of  850nm  to  839nm.  The  blue-shifting  occurs  because  of  the  choice  of 
periodic  structure  of  the  HL  layers.  By  making  the  cavity  2.2X  the  dip  can  be  red- 
shifted  back  to  850nm,  if  desired. 

Figure  39  shows  the  normal  incidence  reflectance  and  phase  calculated  without 
absorptance.  The  most  obvious  difference  is  in  the  minimum  value  of  the  dip.  The 
absorptance  causes  a  perceived  lower  reflectance  on  the  lower  mirror.  This  is  due  to 
less  photons  being  reflected  due  to  absorption.  Therefore  the  dip  is  lower  with 
included  absorptance.  The  linear  phase  region  used  by  Corzine  and  Babic  and 
described  in  Chapter  2  is  clearly  seen  between  830nm  and  850nm  in  both  figures. 


M  m  X  XI2 


20nm  cap 
layer 


Figure  37.  Composition  diagram  for  Structure  2 
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Figure  38.  Calculated  reflectance  and  phase  spectra  for 
Structure  2,  including  absorption,  (k  *  0) 
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Figure  39.  Calculated  reflectance  and  phase  spectra  for 
Structure  2,  not  including  absorption,  (k  =  0) 
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Figure  40  shows  an  example  of  the  linear  periodic  grading  schemes  that  can  be 
used  in  the  HL  half  wave  DBR  layers.  If  Scheme  1  is  used,  the  Fabry-Perot  dip  blue- 
shifts,  unless  the  cavity  is  extended.  If  Scheme  2  is  used  the  Fabry-Perot  dip  is  at 
the  design  wavelength.  Figure  41  shows  the  refractive  index  profile  at  XQ  -  850nm 

for  Structure  2  using  grading  Scheme  1.  The  cavity  region  and  a  single  HL  layer  on 
either  side  of  the  cavity  are  shown.  The  cavity  contains  four  quantum  wells.  The 
stepped  profile  designates  the  sublayers  used  in  the  Mathcad  calculation. 


Figure  40.  Example  periodic  grading  schemes  of  high-low 
index  layers 
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Sublayers  of  RCLED 

Figure  41.  Actual  refractive  index  profile  of  Structure  2  at 
XQ  -  850nm  used  in  reflectance  program 
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Figure  42  shows  the  reflectance  of  the  upper  and  lower  DBRs  seen  from  inside 


the  microcavity  looking  out  for  TE  Polarization. 


Figure  42.  Calculated  reflectance  spectra  for  upper  and 
lower  DBRs  of  Structure  2  at  0°  incidence  angle 

3.4  The  Spontaneous  Emission  Enhancement  Factor:  H 

In  Chapter  2  it  is  stated  that  the  threshold  gain  necessary  for  stimulated 
emission  is  not  achievable  due  to  the  low  reflectance  of  the  upper  DBR.  However, 
the  presence  of  the  microcavity  does  change  the  spontaneous  emission  profile  in  a 
measurable  way.  The  spontaneous  emission  is  enhanced  or  diminished  depending  on 
the  location  of  the  quantum  well  with  respect  to  the  resonant  standing  wave 
anti-nodes  of  the  microcavity.  This  change  is  characterized  by  the  spontaneous 
emission  enhancement  factor  derived  here  using  a  classical  wave  interference  model. 

Presented  here,  for  the  first  time,  is  the  inclusion  of  angular  dependency  in  S. 
The  angular  dependence  represents  the  angular  sensitivity  of  the  spectral 
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enhancement  of  S.  As  the  angle  from  the  normal  increases,  shorter  wavelengths  will 
be  enhanced  and  the  peak  wavelength  of  the  emitted  light  blue-shifts. 

3.4.1  Derivation  of  E.  Figure  43  shows  the  architecture  used  to  derive  S.  The 
microcavity  contains  a  layer  of  emitters  located  at  dt  from  the  upper  DBR.  This  is 

the  QW  layer.  Multiple  QWs  are  assumed  thin  enough  to  be  considered  as  a  single 
layer,  although  in  the  device  the  QW  layer  may  take  up  as  much  as  25%  of  the 
thickness  of  the  cavity.  For  example,  in  a  RCLED  emitting  at  850nm,  a  IX,  cavity  is 
274nm  thick,  using  ncav  «  3.1.  A  RCLED  with  four  QWs,  each  8nm  wide  with  3 
barriers  of  lOnm  represents  a  layer  72nm  thick.  In  this  study,  the  finite  thickness  of 
the  QWs  is  not  considered.  The  QWs  are  assumed  to  be  represented  by  a  layer  of 
infinitesimal  width  in  the  cavity. 

The  DBRs  have  complex  reflectivity  pleJ0'  and  p2eJ0i  which  contain  amplitude 
and  phase  information  necessary  for  calculating  reflected  fields.  The  emitted 
radiation  is  examined  using  a  classical  wave  interference  model. 


Quantum  Well  Layer  = 


Figure  43.  Graphic  representation  of  architecture  for  H 
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For  emission  at  an  angle  (j),  two  rays  may  be  considered  to  propagate.  Ray  1 

■ 

(denoted  by  the  solid  arrow  in  Figure  43)  is  directed  up  towards  the  upper  DBR  and 
is  partially  reflected  and  partially  transmitted  by  that  DBR.  Ray  2  (denoted  by  the 
dashed  arrow  in  Figure  43)  is  directed  down  towards  the  lower  DBR  and  is  partially 
reflected  and  partially  transmitted  by  that  DBR.  Subsequent  reflections  and 
transmissions  through  the  upper  DBR  for  both  rays  are  also  shown  in  Figure  43. 

To  analyze  the  microcavity  the  relative  phase  differences  between  the  ray  and  its 
reflections  ( Supper )  and  the  second  ray  and  its  reflections  {Slower)  must  be  accounted  for. 

Then  the  relative  phase  difference  between  Ray  1  and  Ray  2  is  needed  in  order  to 
interferometrically  sum  their  fields.  The  relative  phase  differences,  {Supper)  and 
{Slower),  are  equal  for  both  rays.  They  will  now  be  referred  to  as  S,  which  is  given  in 

Equation  52  as  a  function  of  X,  the  emitted  wavelength,  and  2,  the  growth  direction. 

The  functional  dependence  on  the  growth  direction,  z,  is  due  to  the  cavity  grading 
shown  in  Figure  37. 

A,  ■rr 

Sx(X,z)  = — n{X,z)d  cos(^)  ( radians )  (52) 

X 

The  modified  E-field  spectrum  ,  Ep  for  Ray  1  in  a  microcavity  is  given  by: 

Ex  (X)  =  r,  Eqw  +  r xp2px  Eqwe-5'  +. .  .+rx (p2px )”  Eqwe~niS'  (Volts  fm)  (53) 

where  t1  is  the  transmissivity  through  the  upper  DBR,  Eqw  is  the  unmodified  E-field 
spectrum,  and  n  denotes  the  number  of  reflections.  By  summing  n  over  infinity  E1 
becomes: 
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£x  (A, z)  —  Tj  (A)  E ^  (A)  i  _  ^ (A)^  ( A) g_, ,, (Xj) 


(Volts/m) 


(54) 


The  E-field  due  to  Ray  2  is  also  modified  and  becomes  E2  given  by: 

E2(A)  =  Txp2Eqw+Txp2pxp2Eqwe"s'+...+Tl(pxp2)n  Eqwe~mS'  (Volts/m)  (55) 


which  reduces  to: 


E2  ( A,z )  =  tx  ( A)p2  (A)  E  (A) 


1 


\-p2(A)px(A)e 


-iS,(X,z) 


(Volts/m) 


(56) 


The  relative  phase  difference,  8  ,  between  Ray  1  and  Ray  2  is  given  by: 


4  71 

8rel(A,z)  =  —n(A,z)d2  cos(0) 


(57) 


where  d2  is  the  distance  from  the  QW  layer  in  RCLED  to  the  lower  DBR.  Using  8rel 
the  modified  E-fields  may  be  summed  interferometrically  to  obtain  the  output 
E-field  measured  at  the  surface  of  the  RCLED,  Eemitted. 


EemittM,z)  =  W)  +  E2We-iS”iaJ) 


(58) 


E emitted  i^’Z)  ~  U 


1  +  p2e~'K,{X'z) 
1  -PM"' 


EqwW 


(59) 


The  output  intensity  is  proportional  to  the  modulus  squared  of  the  output 
E-field,  Eemitted.  The  ratio  of  I  Eemitted  1 2  to  I  Eqw  1 2  represents  the  Spontaneous 
Emission  Enhancement  Factor,  E,  due  to  the  microcavity  and,  for  emission  through 
Rp  is  given  by: 
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2(A  ,z)  = 


J emitted 


Hw\ 


(60) 


which,  after  substitution,  becomes: 


1  + 


2(A,Z>  : 


\p2\2  +  2|  |p2 1  cos^02  (A)  -  ^ n(A,z)(d  -  dx )  cos^ 


1  +  |p,  |2  \p2 12  -  2jpj  || p2 1  cosftf,  (A)  +  02  (A)  -^«(A,z)J cos^ 


(61) 


For  the  case  of  no  absorption  in  the  cavity  or  the  DBRs,  the  transmission,  I  r7 1 2, 
through  the  upper  DBR  is  equal  to  (1  -Rt)  and  Equation  60  becomes: 


2(A,z)  = 


(l  -  R((A))I  1  +  (A) +  2^  (A)  cosj^02  (A)  -  —«(A,z)(<i  -<i,)cos0 
1  +  ^(AWA)  -  2a/R,(A)R2(A)co^,(A)  +  02(X)  -  ~n(X,z)dcos<t> 


(62) 


where  R,  and  R2  are  the  reflectances  of  the  upper  and  lower  DBRs  respectively. 

3.4.2  Applying  E.  The  angle  <|>  in  Equation  60  must  be  referenced  to  the  output 
angle  measured  from  the  surface  of  the  RCLED,  a.  Using  Snell’s  Law  for  refractive 
indices: 

nairsina  =  nqwsm<j>  (63) 


where  n  is  the  refractive  index  of  the  quantum  well  layer  of  the  cavity  and  nai-l . 


qw 


Now  substituting  back  into  Equation  60: 

r  ( 


(l-R,)  l  +  R^V/^cosI  G2-^^(d-d,\nJ  -sm2  a) 
_ v _ v  V] _ 

1  +  RlR2-  2  cos  ex  +  e2  -  4*n(z')  d(nj  -  sin2  a)^2 

v  ^yV) 


2(A) 


J) 


(64) 
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Equation  64  includes  absorption  in  the  DBRs,  but  not  in  the  cavity.  Note  that 
(. Rv  9 j)  and  (R2, 02)  must  be  calculated  from  the  viewpoint  of  the  quantum  well  layer 

inside  the  cavity.  This  requires  reversing  the  order  of  the  matrix  multiplication 
outlined  in  Chapter  2  for  Rr  Figure  42  shows  calculated  upper  and  lower  reflectance 
spectra  for  Structure  2.  There  is  no  dependency  on  carrier  injection  level  included  in 
E,  therefore  the  customary  widening  of  the  emission  spectrum  with  increasing 
injection  seen  in  a  conventional  LED  is  not  expected. 

3.4.3  Examining  the  Effects  of  E.  Figure  44  shows  the  blue  shifting  of  S  with 
increased  angle  from  the  normal.  The  displayed  angle  is  referenced  to  the  emission 
angle  from  the  surface  of  Rl  (the  upper  DBR).  The  blue-shifting  has  a  cosine 
dependence  as  seen  in  Equation  60.  The  wavelength  dependent  reflectance  and 
phase  change  of  Structure  2  are  calculated  using  the  Mathcad  document  in 
Appendix  B.  The  graphic  insert  shows  the  location  of  the  QW  layer  in  the  cavity. 
Figure  45  shows  the  calculated  E  (normalized)  with  calculated  reflectance  spectra  of 
Structure  2  at  0°  and  40°.  The  spontaneous  enhancement  is  predicted  to  occur  at  the 
Fabry-Perot  dip.  Figure  46  shows  the  effect  of  placing  the  QW  layer  at  a  node  of  the 
resonant  mode  located  at  A/4  from  the  lower  DBR.  Emission  through  R}  is  greatly 

inhibited  as  compared  to  Figure  44  where  the  QW  layer  is  placed  at  a  cavity  anti¬ 
node.  Figure  47  shows  the  effect  of  placing  the  QW  layer  at  a  node  of  the  microcavity 
located  A/4  from  the  upper  DBR.  Again  the  spontaneous  emission  is  inhibited  greatly 
compared  to  Figure  44. 
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Figure  44.  E  for  Structure  2  with  quantum  well  located  at 
the  center  of  the  cavity  (TE  polarized) 


Wavelength  (nm) 


Figure  45.  Normalized  E  for  Structure  2  at  0°  and  40°  with 
calculated  reflectance  spectrum(TE  polarized) 
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S  (arb.  units) 


Figure  46.  S  calculated  for  quantum  well  layer  located  XI 4 
away  from  the  center  of  the  microcavity  toward  R2 


Figure  47.  S  calculated  for  quantum  well  layer  located  A/4 
away  from  the  center  of  the  microcavity  toward  R1 
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The  next  step  is  to  modify  the  unenhanced,  calculated  emission  spectrum  of  the 
QW  layer  with  E  in  order  to  calculate  the  expected  output  spectrum  of  Structure  2. 
Figure  48  shows  that  the  spontaneous  emission  blue-shifts  with  angle  and  that  the 
Q  of  the  output,  increases  with  increasing  angle. 

Table  4  shows  the  FWHM,  Apeak  and  Q  (Q  -  Xpeak/FWHM/),  associated  with  Figure 

48.  Table  4  indicates  that  this  model  expects  the  narrowing  of  the  emission 
spectrum  with  increasing  output  angle.  This  narrowing  of  the  FWHM  is  attributed 
to  increased  DBR  reflectance  at  shorter  wavelengths  due  to  greater  differences 
between  the  refractive  indices  of  the  high  and  low  index  quarter  wave  layers. 

It  should  be  noted  that  E  could  also  be  modified  by  other  microcavity  factors  in 
order  to  completely  model  the  coupling  between  the  energy  in  the  field  and  the 
emission  mode  of  the  RCLED.  The  mode  coupling  parameter,  p,  accounts  for  this 
efficiency.  P  is  defined  as  the  “cavity-induced  enhancement  of  the  spontaneous 
emission  coupling  efficiency  into  a  lasing  mode  [27].  ”  The  external  quantum 
efficiency  also  decreases  for  increasing  cavity  Q,  because  fewer  photons  are  emitted 
due  to  increased  reflectance  [33],  These  parameters  are  beyond  the  scope  of  a 
classical  wave  interference  model  and  require  a  rigorous  quantum  electrodynamic 
treatment[8;34].  Therefore  Figure  45  does  not  necessarily  show  the  correct  scaled 
output  with  respect  to  intensity.  This  may  affect  other  results,  such  as  FWHM,  that 
must  be  measured  using  intensity. 
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Figure  48.  Calculated  spontaneous  emission  from  Structure 
2  for  carrier  injection  levels  of  1.5  x  1018  cm  3  and  2.5  x  1018 
cm  3.  T  =  300K 


Table  4.  Calculated  Parameters  of  Figure  48. 


Angle  (deg) 

0 

5 

10 

15 

20 

25 

30 

35 

40 

A^peak 

843.1 

842.9 

842.0 

840.7 

838.9 

836.6 

834.0 

831.1 

828.0 

FWHM(AXvO 

3.832 

3.822 

3.741 

3.673 

3.526 

3.382 

3.14 

2.873 

2.568 

Q 

220 

220 

225 

229 

238 

247 

266 

289 

322 

3.5  Summary 

A  new  classical  wave  interference  model  with  included  angular  dependence  was 
presented.  The  purpose  was  to  examine  the  enhancement  effects  of  a  microcavity  on 
spontaneous  emission.  First,  the  reflectance  spectrum  of  Structure  2  was  calculated 
with  the  Mathcad  document  in  Appendix  B.  Then,  the  spontaneous  emission 
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spectrum  of  Structure  1  was  calculated  using  Appendix  C.  Finally,  the  spontaneous 
emission  enhancement,  S,  due  to  the  presence  of  a  microcavity  was  calculated  for 
Structure  2.  5  enhances  the  spontaneous  emission  of  Structure  2  and  also  blue- 
shifts  the  emission  spectrum,  with  increasing  angle. 

The  choice  of  grading  scheme  used  for  the  DBRs  changes  the  effective  size  of  the 
cavity.  Using  Scheme  1  in  Figure  40  blue-shifts  the  Fabry-Perot  dip  while  Grading 
Scheme  2  puts  the  dip  at  the  calculated  wavelength.  Scheme  1  may  be  favored 
because  of  the  better  field  confinement  for  the  quantum  well  layer  in  the  separate 
confinement  hetero-structure  formed  by  the  cavity. 

Other  E  modification  factors  caused  by  the  micro-cavity  such  as  the  mode 
coupling  efficiency  p,  and  decreasing  external  quantum  efficiency  with  increasing  Q 
are  not  considered  in  this  model. 
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4.  Experimental  Method 


4.1  Laboratory  Configuration 

This  section  contains  the  experimental  setups  used  for  measuring  the  photo¬ 
luminescence  (PL),  electro-luminescence  (EL),  reflectance  and  angular  emission  spectra 
under  EL  of  the  two  cavity  structures.  There  were  two  samples  available  for  study.  Sample 
E236  corresponds  to  Structure  1  and  consists  of  a  microcavity  alone.  Sample  E231 
corresponds  to  Structure  2  and  is  a  complete  RCLED  structure. 


2” 


Figure  49.  Location  of  samples  E236  and  E231  on  2”  GaAs  wafer 

E236  and  E231  both  consisted  of  processed  and  unprocessed  halves  of  a  strip  of  the 
wafer  taken  from  the  center  (Figure  49).  Processed  means  that  metal  contacts  with  circular 
apertures  for  light-emission  were  deposited  on  the  top  and  bottom  of  the  devices  after 
growth  using  lithographic  techniques.  The  apertures  range  in  diameter  from  5pm  to  50pm. 
Unprocessed  means  that  no  fabrication  was  performed  after  growth  (i.e.  bare  epitaxial 
layers). 
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5  pm  10  |xm  20  |xm  30  |xm  40  |xm  50  |xm 


Figure  50.  Device  layout  on  processed  samples  (the  size  of  the 
aperture  is  indicated) 

4.2  Photo-Luminescence  Experimental  Setup 

Photo-luminescence  (PL)  is  caused  by  the  injection  of  energetic  photons.  The  photons 
must  have  energy  Einject  >  Euan,  where  Euan,  is  the  transition  energy  of  the  energy  levels  in 
the  QWs.  The  injected  energy  generates  carriers  that  recombine  at  the  QW  transition 
energies.  For  an  8nm  GaAs  QW  with  Al0 15Ga085As  barrier  layers  this  represents  a  peak 

emission  wavelength  of  ~850nm. 

The  experimental  setup  necessary  to  measure  both  the  emitted  PL  and  reflectance 
spectra  during  PL  is  shown  in  Figure  51. 

The  Ti:Sapphire  laser  is  pumped  by  an  Argon  laser  and  is  wavelength  tunable  from 
~780-920nm.  The  white  light  source  emits  broadband  light  (ultra-violet  -  infrared)  to 
generate  the  reflectance  spectrum.  The  Ti:Sapphire  must  be  tuned  to  a  wavelength  where 
the  RCLED  has  low  reflectance,  and  which  is  less  than  the  design  wavelength,  in  order  to 
pump  the  RCLED  and  fill  the  QWs  with  carriers.  The  white  light  source  is  used  to 
determine  the  reflectance  of  the  substrate  at  the  pumping  location.  To  ensure  that  the 
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Ti:Sapphire  laser  and  the  white  light  source  illuminate  the  RCLED  at  the  same  location,  it 
is  necessary  to  align  them  using  beam-splitters  (Figure  51).  Peak  PL  emission  occurs  at  the 
Fabry-Perot  dip  of  the  RCLED  being  pumped. 


Figure  51.  Experimental  setup  for  photo-luminescence  and 
reflectance  measurements  on  RCLEDs. 


Table  5.  Equipment  Listing  for  Photo-Luminescence. 


Optical  Multichannel 

Manufacturer 

Princeton  Applied  Research 

Analyzer  (OMA) 

Model 

EG&G  1460 

Accuracy 

±3A 

Ti:Sapphire  Laser 

Manufacturer 

Coherent 

Output  Power 

~1  Watt 

Diode  Focusing  Lens 

Manufacturer 

Melles  Griot 

Focal  length 

2.54cm 

Spot  Size 

1.85qm 

4.3  Electro-Luminescence  Experimental  Setup 

Electro-luminescence  (EL)  is  caused  by  the  application  of  an  electric  field.  As  the  field  is 
applied  to  the  RCLED,  injected  carriers  recombine  in  the  QW  layers  and  photons  are 
emitted.  The  wavelength  of  the  emitted  photons  is  governed  by  the  energy  levels  within  the 
QW  of  the  RCLED  microcavity.  The  experimental  setup  necessary  to  measure  both  the 
emitted  EL  and  reflectance  spectra  during  EL  is  shown  in  Figure  52. 


Optical 

Multi- 

Channel 

Analyzer 

(OMA) 

— r~ 


Computer 


with  OMA 
software  and 


J?  "  Monitor 


M  =  Mirror 
BS  =  Beamsplitter 


CCD  Camera 


Probe 


BS2 


|  X  1UUV  /  | 

- V  Qpsc 


f  Ml 


(tscope  objective 
(x,y  stage) 


Device  (x,y  stage) 


Figure  52.  Experimental  setup  for  electro-luminescence  and 
reflectance  measurements  of  RCLEDs 


EL  requires  a  charge -coupled  device  (CCD)  camera  in  order  to  place  the  electrical  probe 
on  the  RCLED  near  the  aperture.  It  was  found  that  imaging  the  devices  directly  onto  the 
CCD  array  without  the  camera  lens  worked  best.  During  this  research  the  electrical  probe 
was  located  at  exactly  the  same  place  on  the  device  each  time  in  order  to  achieve  uniformity 
in  injection  conditions.  The  device  was  also  aligned  in  the  same  orientation  in  the  field  of 
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view  of  the  microscope  objective  to  reduce  errors  due  to  different  sampling  conditions.  The 
microscope  objective  has  magnification  10X  and  numerical  aperture  of  0.25. 

A  waveform  generator  in  continuous  wave  (CW)  mode  was  used  to  energize  the 
electrical  probe  and  cause  electro-luminescence  in  the  RCLED.  The  electrical  path  to 
ground  for  the  injected  electrons  is  through  the  RCLED  and  into  the  x-y  translation  stage 
(grounding  pad)  shown  in  Figure  52.  It  is  essential  that  good  electrical  contact  be  made 
between  the  back  of  the  RCLED  and  the  grounding  pad. 

4.4  Angular  Measurements  Using  Fiber  Optic  Probe 

In  order  to  fulfill  the  main  goal  of  this  thesis,  the  correlation  of  the  theoretical  emission 
spectrum  with  the  experimentally  measured  spectrum,  it  was  necessary  to  build  a  fiber 
optic  probe  capable  of  being  rotated  around  the  sample.  The  silica  fiber  used  was  a  Belden 
Inc.,  graded  index  , multimode  fiber  with  a  50jum  core  and  a  numerical  aperture  (NA)  of 
0.28. 

The  tip  of  the  fiber  probe  was  aligned  to  the  axis  of  rotation  of  the  angular  translation 
device  in  Figure  53.  Therefore  the  tip  of  the  probe  would  not  be  pulled  off  the  RCLED 
during  rotation.  In  practice  it  was  better  to  intentionally  pull  the  fiber  tip  off  the  axis  of 
rotation  by  a  small  amount  (Figure  54),  because  if  it  were  too  long  the  fiber  would  hit  the 
RCLED  during  rotation.  Some  assumptions  were  made  in  order  to  validate  the  use  of  the 
fiber  probe.  First,  that  the  best  way  to  keep  the  fiber  probe  on-sample  (that  is  pointing 
squarely  at  the  sample)  was  to  use  the  intensity  reading  from  the  optical  multi-channel 
analyzer  (OMA)  as  the  metric  of  alignment.  That  is,  for  a  given  angle,  a,  the  fiber  was 
translated  using  its  x-y  stage  until  the  peak  wavelength  intensity  was  maximized, 
according  to  the  OMA.  Second,  that  the  fiber  used  was  too  short  ( ~0.5m )  to  significantly 
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change  the  parameters  being  measured,  either  through  absorption  or  dispersion.  The  fiber 
probe  is  shown  in  Figure  53. 

The  physical  limitations  of  the  fiber  probe  are  as  follows.  The  angle  of  rotation  of  the 
fiber  optic  was  measured  to  an  error  of  ± 1  degree.  The  fiber-optic  probe  was  aligned  to  the 
axis  of  rotation  to  an  error  of  ~-lmm.  Again,  this  was  necessary  to  prevent  the  probe  being 
pulled  off  the  sample  as  it  was  rotated  (Figure  54)and  the  error  was  purposely  made  on  the 
short  side  in  order  to  keep  the  probe  tip  on  a  constant  circular  trajectory  as  it  was 
translated  around  the  sample  This  kept  lateral  translations  with  the  x-y  stage  to  a 
minimum. 

The  multi-mode  nature  of  the  fiber  caused  broadening  in  the  measured  samples.  This 
broadening  is  due  to  two  factors.  First,  the  numerical  aperture  (NA)  caused  emission  from 
non-desired  angles  (not  the  angle  being  measured)  to  enter  during  a  measurement  (Figure 
54).  Secondly,  due  to  the  multi-mode  nature  of  the  fiber,  the  emission  from  the  non-desired 
angles  propagated  in  the  fiber  and  broadened  the  measurement.  The  broadening  is  less  for 
lesser  angles  (~0°)  because  of  the  symmetry  of  the  emission  at  those  angles.  As  the 
measured  angle  increases  the  asymmetrical  field  at  higher  angles  causes  more  broadening 
in  the  measurement.  Broadening  of  the  emission  spectrum  leads  to  errors  in  measured  full 
width  at  half  maximum  (FWHM).  Fiber  broadening  would  be  avoided  with  a  single-mode 
fiber,  however  the  equipment  was  not  available  to  align  a  single-mode  fiber  with  ~  7gm 
core. 
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Figure  53.  Fiber-optic  probe  with  angular  translation  used  in 
experiments. 
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Figure  54.  Fiber  being  rotated  around  sample  RCLED,  the  dashed 
lines  represent  the  numerical  aperture  of  the  fiber 
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4.5  Aligning  The  Fiber  Optic  Probe 

The  best  way  to  align  the  fiber  probe  was  to  first  electrically  pump  the  RCLED  using 
the  microscope  objective  to  image  the  electrical  probe  and  the  RCLED  on  the  CCD  camera. 
Then  the  microscope  objective  was  removed  and  the  fiber  probe  put  in  its  place.  The  output 
end  of  the  fiber  probe  is  placed  against  the  CCD  camera.  When  the  fiber  probe  is  aligned  a 
circular  pattern  appears  indicating  a  photon  flux  emerging  from  the  fiber-optic.  Next  the 
end  of  the  fiber  probe  is  placed  against  the  entrance  slit  of  the  OMA  and  adjusted  until  an 
emission  pattern  is  visible  on  the  computer  monitor.  The  fiber  is  now  vertically  aligned  with 
the  sample  and  ready  for  angular  measurements.  The  probe  worked  extremely  well. 
Although  it  was  always  difficult  to  align  initially,  once  aligned  it  was  relatively  easy  stay 
centered  on  the  sample. 
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5.  Experimental  Results 


5.1  Introduction 

For  the  purposes  of  presenting  the  experimental  data,  sample  E236  refers  to 
Structure  1  and  sample  E231  refers  to  Structure  2.  Both  designs  have  been  discussed  in 
Chapter  3.  As  shown  in  Figure  55  both  E236  and  E231  have  processed  and  unprocessed 
sections.  “Processing”  means  depositing  metal  contacts  on  the  upper  and  lower  wafer 
surfaces  in  order  to  electrically  inject  carriers  into  the  sample.  Photo-luminescence  does  not 
require  these  extra  manufacturing  steps.  The  wafer  had  been  split  into  three  sections  for 
each  sample.  One  section  was  unprocessed  and  was  used  for  PL.  Two  sections  were 
processed  and  used  for  EL.  Sections  ii  and  iii  had  very  similar  properties  because  the 
sample  was  rotated  during  growth.  Section  iii  was  not  used  during  the  experiment. 


5  cm/2  inches 

-4 - ► 

i  ii  iii 


As  grown  Fabricated 

Figure  55.  Diagram  of  E236  and  E231  orientation  used  in  PL  (left) 
and  EL  (right) 

This  chapter  presents  the  data  measured  using  the  methods  outlined  in  Chapter  4  and 
then  compares  these  measured  results  to  the  theoretical  results  of  Chapter  3.  Good 
agreement  is  found  for  predicting  blue-shift  behavior  with  angle.  The  FWHM  is  not  tracked 
very  well  due  to  the  limitations  of  the  model  with  regard  to  changing  absorptance  under 
current  injection,  lack  of  second-order  microcavity  effects,  and  diffraction  effects  due  to  the 
device  aperture. 
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5.2  Device  Parameters 


These  devices  were  built  using  a  low-pressure  MOVPE  system  with  quartz  lamp 
heating.  The  temperatures  inside  the  growth  chamber  can  reach  700-800°C  depending  on 
the  layer  being  grown.  The  layers  of  E231  and  E236  are  grown  on  top  of  an  ~500/im  thick, 
(100)  oriented,  (n+)  GaAs  substrate.  A  buffer  layer  of  (n+)  GaAs  was  grown  on  the 
substrate  to  smooth  the  surface  for  further  growth.  E231’s  lower  DBR  (Si  doped  ~ 
2xl018cnr3)  has  a  HL  index  X/2  period  oi~12nm  linear  grading  A^Ga^As:  0.15<x<1.0, 

~58.6nm  AlAs,  ~12nm  linear  grading  AlxGa1.xAs:  1.0>x>0.15,  ~49.8nm  Al015Ga0  85As.  There 
are  38.5  HL  periods  in  the  lower  DBR.  E231’s  upper  DBR  (C  doped  ~  2xl018cnr3)  has  a  LH 
index  X/2  period  of  ~58.6nm  AlAs,  ~12nm  linear  grading  AljcGa1.xAs:  1.0>x>0.15,  ~49.8nm 
Al0 15Ga085As  and  ~12nm  linear  grading  AlxGa1.xAs:  0.15<x<1.0.  There  are  6  periods  in  the 

upper  DBR.  The  last  X/4  layer  contains  ~20nm  p+  GaAs  to  help  form  an  ohmic  contact.  The 
E231  DBR  grading  scheme  corresponds  to  Scheme  1  (Chapter  3)  and  the  Fabry-Perot  dip  is 
expected  to  blue  shift  from  the  design  wavelength  of  850nm.  E236  has  ~735nm 
Al0 15Ga0  85As,  ~12nm  linear  grading  AlxGal  xAs:  0.15<x<1.0,  ~58.6nm  AlAs  beneath  its 

cavity  and  ~58.6nm  AlAs,  ~12nm  linear  grading  AlxGal  xAs:  1.0>x>0.15,  ~715nm 
Al0  i5Ga0  85As,  ~20nm  p+  GaAs  above  the  cavity.  The  cavity  for  both  E236  and  E231  is 
~99.8nm  AlxGa1.xAs:  1.0>x>0.15,  four  GaAs  8nm  quantum  wells  with  lOnm  Al015Ga085As 
barriers,  ~99.8nm  AlxGa1.xAs:  0.15<x<1.0.  The  cavity  is  IX.  thick  at  X  =  850nm.  This 
information  is  abbreviated  in  Table  6. 
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Table  6.  Growth  Table  for  E231  and  E236. 


E236  (Structure  1) 

E231  (Structure  2) 

Ohmic  Contact 

200A  GaAs  Cap  p+doped  (C) 

200A  GaAs  Ca  p+doped  (C) 

Upper  DBR  (HL) 

7150A  Buffer  Region  p-doped  (C) 

Al0 15Ga0  85As/A1As  (linearly 
graded)  p-doped  (C) 

Cavity 

IX  Al0 15Ga0  85As/A1As  (linearly 
graded) 

IX  Al0 15Ga0  85As/A1As  (linearly 
graded  over  998  A) 

Quantum  Wells 

Four  8nm  GaAs  QWs  with  lOnm 
Al015Ga085As  Barriers 

Four  8nm  GaAs  QWs  with  lOnm 

A1q  15Ga0  85As  Barriers 

Lower  DBR  (HL) 

7350A  Buffer  Region  n-doped  (Si) 

Al0 15Ga0  85As/A1As  (linearly 
graded)  n-doped  (Si) 

Substrate 

GaAs  n+  doped  (Si) 

GaAs  n+  doped  (Si) 

5.3  Experimental  Data  Presentation 

Figure  56  shows  the  measured  blue-shifting  of  the  emission  spectrum  with  increasing 
angle  of  E231  (Structure  2)  under  electro-luminescence.  This  was  expected  from  the 
calculations  in  Chapter  3.  The  spectrum  broadening,  instead  of  narrowing,  with  increasing 
angle  is  attributed  to  diffraction  from  the  50{im  aperture,  sampling  with  a  multi-mode  fiber, 
and  increased  absorptance  with  injected  current.  Diffraction  causes  energy  from  smaller 
emission  angles  to  be  spread  into  larger  measured  angles.  Other  apertures  were  not  used 
for  measuring  because  of  difficulties  in  aligning  the  fiber  probe. 

Figure  57  shows  the  calculated  emission  from  E231  (Structure  2)  using  the  measured 
emission  spectrum  of  E236  (Structure  1).  The  peak  wavelengths  and  FWHMs  of  Figure  57 
are  given  in  Table  4. 
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Figure  56.  Measured  electro-luminescence  spectrum  of  E231 
(Structure  2)  versus  angular  deviation  from  normal 


820  825  830  835  840  845  850  855  860 


Wavelength  (nm) 

Figure  57.  Normalized  calculated  E231  (Structure  2)  emission 
(solid)  using  measured  E236(Structure  1)  emission. 
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Table  7  contains  the  relevant  parameters  from  Figure  56.  These  parameters  are  peak 
wavelength,  FWHM,  and  Q  versus  measured  angle.  Figure  58  shows  the  measured 
reflectance  from  the  unprocessed  E231  sample.  The  measurements  were  taken  at  the 
corresponding  locations  in  Figure  55.  The  Fabry-Perot  dip  does  not  vary  by  more  than  1% 
between  devices  4  and  6  showing  high  uniformity  during  growth.  Device  2  is  ~3%  from  the 
other  devices,  probably  because  it  was  located  so  close  to  an  edge,  which  can  cause 
non-uniformity.  Figure  59  compares  the  calculated  and  measured  reflectances  for  E231 
(Structure  2).  Figure  59  shows  that  the  calculated  reflectance  does  not  include  enough 
absorptance  in  the  refractive  index  information.  This  is  evident  because  of  the  lower 
reflectances  of  the  measured  data.  The  design  wavelength  is  850nm,  but  the  MATH  CAD 
reflectance  program  calculates  a  shifted  dip  at  841nm  due  to  grading  the  DBRs  using 
Scheme  1.  Figure  60  shows  the  emission  spectrum  of  E236  under  electro-luminescence.  The 
emission  clearly  has  no  angular  dependence  and  does  not  change  for  various  fiber  sampling 
angles.  It  has  the  characteristics  of  a  lambertian  emitter. 


Table  7.  Measured  Parameters  for  Figure  56. 


Angle 

0 

5 

10 

15 

20 

25 

30 

35 

40 

Peak  X 

847.62 

847.32 

846.47 

845.68 

844.65 

842.17 

839.5 

837.08 

834.47 

FWHM 

4.52 

m 

5.07 

5.48 

5.9 

6.53 

6.6 

6.3 

6.84 

Q 

187.53 

180.28 

166.96 

154.32 

143.16 

128.97 

127.20 

132.87 

122.00 
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Reflectance 


Device  2 


740  760  780  800  820  840  860 

Wavelength  (nm) 


880  900  920 


Figure  58.  Measured  reflectance  from  various  locations  on 
unprocessed  E231  sample  (Structure  2) 


Wavelength  (nm) 


Figure  59.  Measured  (solid)  and  calculated  (dashed)  reflectance 
from  device  E231  at  location  3  (unprocessed)  in  Figure  55 
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The  emission  from  a  QW  layer  such  as  E236  (Structure  1)  should  be  narrower  than 
that  of  a  bulk  emitter.  The  spectral  width  of  the  emission  in  Figure  61,  using  the  full-width 
at  half  maximum  (FWHM)  criteria,  is  ~25nm.  The  approximate  predicted  wavelength  for  a 
bulk  spontaneous  emitter  at  300K  with  Xpeak  ~840nm  is  27nm  using  AA  «  1  A5Z2peakkBT , 

where  AX.  and  Apeak  are  in  microns  and  kBT  is  in  eV.  The  QW  emission  is  slightly  narrower 
than  the  calculated  bulk  emission.  Table  2  shows  the  parameters  of  peak  wavelength, 
FWHM,  and  Q  for  Figure  61.  The  FWHM  is  seen  to  increase  fairly  substantially  with 
increased  injection  current.  The  increase  in  FWHM  is  expected  because  of  the  increasing 
separation  of  the  quasi-fermi  levels  with  increased  carrier  injection. 


5-7 


3000 


Figure  61.  Measured  electro-luminescence  spectrum  from  sample 
E236  (Structure  1)  for  various  injection  currents 


The  lack  of  additional  transitions  in  the  spectrum  enforces  the  calculated  finding  of  only 
one  electron/heavy-hole  transition  at  X=850nm  in  the  well.  The  shifting  of  the  transition 
wavelength  for  E236  (Structure  1)  from  the  calculated  850nm  may  be  due  to  thickness 
variations  in  the  quantum  well.  Because  the  transition  has  blue-shifted  from  850nm  this 
indicates  that  the  quantum  well  is  narrower  than  designed,  thus  causing  the  allowed, 
quantized,  energy  levels  to  increase  their  separation.  Figure  63  shows  the  emission 
spectrum  from  a  device  on  E231  plotted  for  various  injection  currents.  Notice  that  the 
FWHM  of  ~5nm  does  not  change  for  increasing  carrier  injection.  This  is  due  to  the 
microcavity  effects  of  the  RCLED.  Although  direct  intensity  comparisons  between  Figures 
61  and  62  is  not  possible  because  no  exact  intensity  readings  were  taken,  a  qualitative 
examination  of  the  numbers  reveals  that  the  spontaneous  emission  has  been  greatly 
enhanced.  Much  narrower  FWHM  and  higher  output  intensity  are  the  main  advantages  of 
RCLEDs  over  conventional  LEDs. 
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Intensity  (arb.  units) 


Table  8.  Measured  Parameters  for  Figure  61. 


Current  (mA) 

20 

30 

40 

50 

60 

Peak  X 

840.8 

841.2 

841.2 

841.4 

841.2 

FWHM 

20.1 

23.0 

25.0 

26.5 

28.0 

Q 

41.8 

36.6 

33.6 

31.8 

30.0 

840  841  842  843  844  845  846  847  848  849  850  851  852  853  854  855 

Wavelength  (nm) 

Figure  62.  Measured  electro-luminescence  of  sample  E231 
(Structure  2)  for  various  injection  currents 


Table  9.  Measured  Parameters  for  Figure  62. 


Current  (mA) 

13 

15 

17 

19 

21 

Peak  X 

847.3 

847.5 

847.6 

847.6 

847.6 

FWHM 

5.2 

5.0 

5.0 

5.0 

4.8 

Q 

162.9 

169.5 

169.5 

169.5 

176.6 

Figure  63  shows  the  emission  and  reflectance  spectra  for  a  device  on  sample  E231. 
Figure  63  clearly  shows  that  the  Fabry-Perot  dip  does  indeed  define  the  location  the 
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emitted  radiation.  Figures  64  and  65  show  the  emitted  spectrum  from  samples  E236  and 
E231  under  photo-luminescence  from  a  Ti:Sapphire  laser.  The  AOM  scale  is  not  power- 
calibrated  but  a  relative  comparison  is  possible.  The  emission  wavelength  of  the 
Ti:Sapphire  laser  is  800nm.  The  same  characteristics  of  FWHM  and  Q  as  for  EL  are  visible, 
although  Figure  64  has  a  lot  of  noise  that  may  be  due  to  photo-luminescence  of  the  bulk 
GaAs  substrate  and  AlGaAs  buffer  by  the  Ti:Sapphire  laser. 


w 

C 

<L> 


Figure  63.  Measured  electro-luminescence  and  reflectance  spectra 
from  same  location  on  sample  E231  (Structure  2) 
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Figure  64.  Photo-luminescence  of  sample  E236  (Structure  1)  at 
various  input  optical  power  densities 
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850 


Figure  65.  Photo-luminescence  from  sample  E231  (Structure  2) 


5-11 


Reflectance 


5.4  Comparison  of  Calculated  and  Experimental  Results 

The  data  in  the  tables  above  is  presented  in  graphical  form  in  direct  comparison  with 
the  calculated  results.  Injection  levels  cannot  be  compared  directly  because  they  were  not 
explicitly  addressed  in  the  model.  Figure  66  shows  the  calculated  Q  and  FWHM  for  the 
plots  of  Figure  61.  The  RCLED  model  predicts  a  non-varying  Q  and  FWHM  with  increased 
injection  current.  This  is  because  the  model  depends  only  on  the  reflectances  of  the  upper 
and  lower  DBRs  of  the  RCLED.  The  Q  has  a  linear  coefficient  of  R2=0.8366.  The  FWHM 
has  a  linear  coefficient  of  Rl=0.823  with  error  bars  denoting  ±3A.  This  error  is  based  upon 
the  accuracy  of  the  OMA  and  the  subtraction  of  the  upper  and  lower  half  maximum 
wavelengths  needed  to  get  the  FWHM.  Figure  66  shows  the  calculated  and  measured 
FHWM  and  Q  for  sample  E231  (Structure  2).  The  parameters  are  plotted  against  the  angle 
<(>,  which  is  the  angle  of  the  fiber-optic  probe  during  measurement.  Figure  67  shows  the 
calculated  and  measured  peak  X  values  for  E231  (Structure  2).  The  values  are  always 
within  1.0%  of  each  other  and  have  a  joint  correlation  coefficient  of  0.998.  The  joint 
correlation  coefficient  is  calculated  using: 

~2l(Xi  -  PxXyt  -  Py) 

p  =  -  -\<p  <\  (unitless)  (65) 

where  u  and  u  are  the  means  of  the  data,  and  cr  and  cr  are  the  respective  standard 

•  x  r  y  x  y 

deviations.  Note  that  the  peak  wavelength  difference  is  due  to  measuring  the  emission  at  a 
location  where  the  DBR  thickness’  and  cavity  thickness  were  not  exactly  as  calculated.  The 
calculated  data  can  easily  be  shifted  to  match  the  experimental  data  by  changing  the 
design  wavelength  used  in  the  emission  model. 
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Figure  68  shows  the  calculated  and  measured  FWHM  for  Structure  2.  RCLED.  The 
trendline  (dotted)  shows  the  coefficient  of  linearity  ,R2,  of  the  experimental  FWHM.  There 
is  a  substantial  disagreement  in  the  change  of  FWHM  vs.  <(>.  The  two  data  sets  have  a  joint 
correlation  coefficient  of  -0.867.  The  broadening  of  the  experimental  FWHM  with 
increasing  angle  may  be  due  to  a  change  of  complex  refractive  index  under  current 
injection.  Increased  absorption,  for  instance  due  to  prolonged  current  injection  during 
measurements,  will  broaden  the  Fabry-Perot  dip  and  cause  the  FWHM  to  decrease. 
Broadening  also  occurs  due  to  diffraction  effects  at  the  device  aperture  and  using  a  multi- 
mode  fiber  for  taking  measurements.  Both  of  these  cause  energy  at  longer  wavelengths 
from  lower  emission  angles  to  be  measured  at  a  higher  angle.  The  diffraction  affects  at  the 
aperture  are  unavoidable.  By  taking  all  the  measurements  at  the  largest  aperture, 
broadening  due  to  diffraction  was  made  as  small  as  possible. 


Figure  68.  Calculated  and  experimental  FWHM  for  E231 
(Structure  2) 
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5.5  Results  Summary 

The  experimental  and  calculated  values  of  FHWM  and  Xpeak  were  examined  for  E231 
(Structure  2).  Good  results  were  obtained  for  tracking  and  estimating  A,peak  with  high 
accuracy,  given  the  initial  location  of  the  Fabry-Perot  dip.  However  the  FWHM  estimator 
shows  high  negative  correlation.  This  may  be  due  to  inaccurate  absorptance  data  under 
current  injection,  diffraction  affects  at  the  device  aperture,  non-inclusion  of  second  order 
microcavity  effects  such  as  mode  coupling  and  loss  of  external  quantum  efficiency,  or  a 
mixture  of  all  four.  All  events  would  serve  to  broaden  the  output.  The  model  does  show  good 
accuracy  in  calculating  the  prospective  range  of  the  output  FHWM,  and  may  prove  a  useful 
tool  in  this  regard. 
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6.  Conclusion 


The  goal  of  this  thesis  was  to  develop  a  to  simulate  the  emission  spectra  of  the  RCLED 
and  its  variation  with  angle.  Achieving  this  goal  required  modeling  the  quantum  well 
energy  levels  and  emission  characteristics,  calculating  gain  and  spontaneous  emission 
spectra  for  the  quantum  well,  calculating  DBR  reflectances,  calculating  electric  field 
intensity,  and  using  a  new  classical  wave  interference  model  to  derive  an  analytic 
expression  for  the  spontaneous  emission  enhancement  factor  (H)  of  a  microcavity  device. 
Two  devices  were  designed,  grown  and  fabricated  in  order  to  measure  the  parameters 
described  above.  The  first  consisted  of  a  IX  cavity  with  four  quantum  wells.  The  second 
consisted  of  the  same  IX  cavity  with  four  quantum  wells,  but  with  DBRs  used  to  form  a 
true  RCLED  structure.  The  emission  from  the  second  device  was  significantly  narrower  in 
bandwidth  and  of  higher  intensity  than  that  of  the  first  device.  The  spectral  emission  of  the 
second  device  was  compared  to  that  of  the  model. 

The  results  achieved  during  this  work  are  new  and  demonstrate  the  emission  pattern  of 
a  RCLED  is  not  lambertian.  A  blue-shift  occurs  in  the  emission  spectrum  with  increasing 
angle  from  the  normal.  The  model  is  extremely  accurate  in  tracking  peak  emission 
wavelength  variation  with  angle.  The  model  was  less  accurate  in  calculating  the  FWHM  of 
the  emitted  radiation,  but  factors  not  considered  in  the  model  come  into  play  for  this 
parameter.  Diffraction  effects  from  the  emission  aperture,  and  the  use  of  multi-mode  fiber, 
cause  broadening  because  longer  wavelengths  from  lower  emitted  angles  are  measured  by 
the  probe  at  a  greater  angle.  Extra  absorption  in  the  device  due  to  heating  under  current 
load  would  also  cause  a  broadening  of  the  emission  spectrum  at  all  angles.  An  experimental 
technique  developed  during  the  course  of  this  work  was  to  use  a  fiber  optic  cable  on  a 
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rotational  stage  to  measure  the  emitted  spectrum  versus  angle  of  an  RCLED  under  electro¬ 
luminescence. 

The  main  contribution  of  this  work  is  a  new  model  for  calculating  the  emission  spectra 
of  the  RCLED  with  angular  dependence.  This  model  is  applicable  to  any  RCLED  structure, 
regardless  of  emission  wavelength  or  DBR  design.  The  Mathcad  document  in  Appendix  B 
can  calculate  the  reflectances  of  any  RCLED  DBR,  regardless  of  compositional  grading 
scheme  or  material  composition,  that  is  required  for  the  model. 

Future  work  in  this  area  should  include  the  following: 

1.  Measure  emitted  spectrum  from  range  of  aperture  sizes  and  compare  using  FWHM' 
as  a  metric.  Diffraction  effects  should  be  greater  for  smaller  apertures. 

2.  Use  a  single-mode  fiber  to  make  angular  measurements.  This  will  eliminate 
broadening  due  to  the  fiber  probe. 

3.  Further  examine  differences  in  DBR  grading  schemes  presented  in  Chapter  3  with 
regard  to  exact  placement  of  the  Fabry-Perot  dip.  This  seems  to  be  a  new  area  for 
research  that  is  presently  untapped. 

4.  Examine  other  microcavity  effects,  such  as  the  mode  coupling  efficiency  factor  /?,  that 
are  not  discussed  in  this  thesis.  /?  modifies  5  in  a  nonlinear  way  and  prevents  all  the 
available  energy  from  being  coupled  into  a  microcavity  emission  mode. 

5.  Apply  model  to  RCLEDs  in  other  wavelength  ranges,  such  as  560-690nm,  using 
(A1  Ga,  )  In,  P  cavities  with  AlGalnP  and/or  AlGaAs  DBRs.  Also  to  the  new  ZnSe 

v  x  1-x'y  L-y 

and  AlGalnN  blue-emitting  RCLEDs 

RCLEDs  are  a  new  and  exciting  addition  to  Photonics,  and  it  was  exciting  to  work  on 
the  cutting  edge  of  LED  technology. 


6-2 


Appendix  A:  Reflectance  Calculations 


A.1  Reflectance  Calculations 

Macleod’s  matrix  method  relating  the  input  and  output  electric  fields  is  used  to 
calculate  the  reflectance  spectrum  for  arbitrary  layers  of  thin  films  [10;26] .  By  relating  the 
electric  and  magnetic  fields  at  two  boundaries  in  a  matrix  form,  it  is  possible  to  calculate 
the  reflectance  spectra  of  complicated  structures. 

Boundary  1 

Boundary  2 

Boundary  3 
Boundary  4 

Boundary  5 

Substrate 

Figure  69.  Plane  wave  incident  on  a  set  of  thin  films 


The  electric  (E)  and  magnetic  (H)  fields  at  boundary  1  may  be  related  to  the  electric  and 
magnetic  fields  at  boundary  2  by: 
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where  Si  is  given  by  Equation  5  for  non-absorbing  media  and  by  Equation  14  for 
absorbing  media.  The  optical  admittance  777  is  measured  in  Siemens  and  is  polarization 
dependent. 
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is  the  characteristic  optical  admittance  of  free  space  and  has  a  value  of  2.6544  x  103 
Siemens.  The  parameter  0J  is  the  angle  between  the  wavefront  and  the  normal  to  the  film 
in  the  first  layer.  For  absorbing  media  the  expressions  for  optical  admittance  change  to 


r/,  =VF (a/,2  -k]-hI  sin2  #0  -  2m,  k]  )1/2  ( Siemens )  (s-polarized  light)  (69) 

7,  =  ¥  — - ^ — - -TjT  ( Siemens )  (p-polarized  light)  (70) 

(«2  -  k]  -  n\  sin2  0O  -  2m, /r,  )1/2 


Using  this  matrix  form  the  fields  in  the  incident  layer  can  be  related  to  the  fields  at  the 
substrate  interface  by: 
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The  reflectance  of  a  series  of  thin  films  is  given  by: 
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3  In  all  equations  an  asterisk  (*)  denotes  the  complex  conjugate. 
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where  Y  is  the  input  optical  admittance  of  the  thin  films  and  is  equal  to  H1/E1.  The 
parameter  Y  can  be  manipulated  into: 

Y-B/C  (Siemens)  (73) 

where  B  and  C  are  obtained  from: 

vC 

The  substrate  material  is  taken  to  be  GaAs  throughout  this  thesis.  The  transmittance 
(T)  and  absorptance  (A)  are  also  calculated  by  Macleod  and  are  readily  accessible  from  the 
same  calculation. 
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APPENDIX  B:  Mathcad  Reflectance  Calculation  Document 


This  document  calculates  reflectance  values  over  a  range  of  wavelengths  for  grading  Scheme  t  structures. 
Output  includes  reflectance  of  the  entire  structure,  reflectance  of  the  upper  and  lower  DBRs  from  a  viewpoint 
inside  the  cavity,  and  the  electric  filed  intensity  at  any  point  in  the  structure  for  any  wavelength  in  the  user 
defined  range.  The  document  uses  Macleod  notation  throughout  and  is  based  on  work  by  Macleod  and  Yeh. 


Constants  used  in  Program 


Hm^m  10  6  nm=m-10'9 


.  ,.-7  newton 

10 - - 

amp 


A^mlO"10  ^=2.6544- 10"3- siemens 

.  s  8.854 1 878 1 7- 1  O' 12-  s  299792458  — 

£  °  me  sec 


eV=joule- 1.60217733  10" 19  h=6.6260755- 10" 34  joule-sec  nair  :=1 


Equations  to  change  wavelength  to  energy  and  vice  versa 

1  74- 1 03 

X(Energy)  :  =  — - -nm  E( wavelength)  :  = 

Energy 

\  eV  j 


1.24- 103 
wavelength 


■eV 


Program  User  Input  Starts  Here 


The  bragg  wavelength:  desired  wavelength  of  operation. 


data„ 


nm 


)u/4  := — -  Jii/4  212.5  *nm  Free  space  quarter  wavelength. 

E(X  opj  =  1.459*eV  The  energy  equivalent  of  the  desired  wavelength. 

Refractive  Index  Calculation  starts  Here.  All  files  use  Aspnes  and  Adachi  AIGaAs  information 

known  :=  (0  0  0  0.099  0.099  .198  .198  .315  .315  .419  .419  .491  .491  .590  .590  .7  .7  .804  .804  1.0  1.0  )T 
A  : = RE ADPRN( ALGAAS )  n  :=  1 ..  cols(A)  -  1  cols(known)  =  1  cols(A)  =21 
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vx 


:  =  A<0>  vy^  :=A<*1>  vs<n>  .  =  lspline(vx,  vy^5-)  f(x,n)  :=interp(vs<i'>,vxJvy<i’>,x 


N  subiayer  :=  19  The  number  of  sub-layers  per  layer 

N  cavity  =51  The  number  of  sublayers  in  the  cavity  stack 

layers  penocj  -  2  •  N  sublayer  layers  penod  “38 


datEj  -N  sublayer 
data2  =N  upper 
data3  -N  cavity 
data4  :  =  N  lower 


z  period  :=1-  layers period  subs  =  layers  period  +  1  air  :  =  0  subs -39 

The  A^Ga^^As  refractive  index  information  is  imported  from  ALGAAS_N.PRN,  and  was  generated  in 
another  file.  The  AlxGal-xAs  index  information  is  taken  from  a  paper  by  Adachi.  Journal  of  Applied 
Physics,  VOL.  58,  No  3, 1  August  1985. 
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nAlGaAs(x>X)  :  =  real 


I  E(X) 


imag 


E(X) 
eV  X; 


Calculated  complex  refractive  Index 


\  eV  x/ 

Refractive  Index  of  low  index  layer 

X  :  =  X lowmdex  +  X highindex  Compositional  average  used  in  calculations  x  =0.575 
avg  2  & 

I x  lowindex  “  x  highindex 


difif 


Compositional  difference  used  in  calculations 


Here  I  calculate  the  actual  high  and  low  indices  for  reference  purposes 

n low  :  =  n  AlGaAs(x  lowindex op)  n  low  =3  001  "  3'379‘ 10  1  data5 :=  n  low 
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nlow  =  3.001  -  3.379-10  10i 


nhigh  nAlGaAs(vxhighindex’AopJ  nhigh 


■  O  ^=3.502-0.0031 


data6:=nhigh  “high  =  3-502 -0.0031 


My  new  approach  is  to  use  a  different  program  for  each  type  of  grading  scheme.  I  do  not  like 
unnecessary  repetitiveness,  however  this  will  simplify  computation  and  speed  up  development.  This 
program  deals  exclusively  with  RAMP  index  profiles.  See  fig  1 . 


tri(x)  :=if(  |x|  <1,1-  |x|  ,0)  Define  tri  function  to  use  in  calculation  of  period  structure,  pi  and  p2 

define  the  center  points  of  the  tri  functions. 


pi  ceH  (N  sublayer)  P*  ^ 

Composition  of  period. 

z  -  pi 


P2  •  2-N  sublayer 


tri(x)  -  tri(6  x)' 


1 


rhomb(x)  ;=- 


x(z)  :  = 


xcompz 


rhomb 


period 


^  sublayer 
"  x^z  period 


x  highindex  x  lowindex) 


T  x  lowindex 


j  -39.44  xcomp.  xlowindex  layers  penod  layers  period  ^  ^ 
2  period  :  =  1  -  44  ncomPz  period  :  =  n  AlGaAs(xcomPz  period’ 1  op) 


subs  :  =  layers  penod+  1 
^  sublayer  -^sublayer  +■  ^ 


z  period > z  period  +  ^y^8  period 


Now  I  will  create  the  cavity  profile  using  the  same  finite  element  method.  I  will  assign  each  of  the  wells  to  a 
particular  sublayer  and  put  barrier  material  between  them.  I  will  then  grade  the  cavity  and  ensure  that  the 
entire  cavity  comes  out  to  measure  the  required  number  of  half-wavelengths. 
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num  :  =  2.0  number  of  half  wavelengths  across  cavity,  may  be  multiple  of  a  half. 

p3  defines  the  center  of  the  cavity  sublayers.  I  then  use  this  to  center  the  rhomb  function  in  the 
cavity.  The  rhomb  function  is  used  to  grade  the  filler  region  of  the  cavity.  The  quantum  wells  and 
barrier  region  are  contained  in  the  fiat  region  on  top. 


p3  :  =  ceil 


N 


cavity 


p3=26  rect(x)  :=if  |x|  ,  1 ,0 


rhomb(x)  :  = 


zcav'  1  ■ 14  cavity 


tri(x)  -  tri(7  x) 

6 
7 


1 


num  qw  ' =  4 
xqw  :  =  0.0001  xbar  •  1 5 


num  barrier  ’  3 


qw:=0..numqw- 1  bar  :=0..  num  barrier-  1 


xcav(z)  =  rhomb  I 


z  -  p3 


N 


cavity' 


]_ 

2/ 


xhighindex  x  lowindex)  x  lowindex  xcavityz  cav  xcav(zcav 


xcavity23^2(|w  xqw  xhighindex  015 


ncavlt>z cav  :  =  n  AlGaAs  xcavityZ cav.^ 


OP 


n  qw  n  AlGaAs 


(xqw’^op) 


L  qw  ;  =  8nm  Barrier  qw  :=  10nm 


The  next  equation  accounts  for  linearly  grading  the  cavity  and  calculates  the  sublayer  thicknesses  at  the 
bragg  wavelength 


cav 


nUm  2~  ~  nUm  clw^  qw^e(n  qw)  num  barrier’ 13  arner  qw  (n  high 
0  ZRe(ncavity)  -  num  qW  Re)n  -  num  barrier  (n  high) 


cav  step  =  44001  *A 


:=L , 


sublayer cavltyz^:  =  cav step  sublayer  =Bamerqw  sublayer  cavity^.^  •  -qw 


sublayer  cavibr  ncavityZcav  =  1  -0.001  i  *1-^  num  A  check  to  ensure  that  the  cavity  measures 
Zcav  cav  \  2  /  up  to  the  proper  number  of  half  wavelengths. 


J  cav 
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User  input  to  decide  how  wide  the  plateaus  of  the  indices  should  be. 
From  all  of  my  design  drawings  it  is  important  to  be  able  to  specify 
this  number  and  then  calculate  the  inter-plateau  distances. 

dkt  j^gh  =  48.674 ’nm  distjow  =  58.799  *nm 


X 

~~  ~  Re(n  high) ' dlst  high  “  Re(n  low) ' dist  low 
StCP  ^aYers  period 

X!  Re(ncomp)n  -  7  Re(nlow)  -  7  Re(nhlgh) 
n=  1 


The  step  size  is  determined  by 
subtracting  the  optical  distance 
(refractive  index  times  physical 
distance)  of  the  plateaus  from  the 
half  wave  layer.  The  remainder  of 
course  is  what  is  left  for  grading. 


a  step  =  801729  *A 


Now  1  make  the  assignations  of  sublayer  size  based  on  the  sizes  disthigh  and  dist,ow. 

sublayerz  p0n0(|  =  a  step 

p2  =  38 

j  =0..  6 

dlst  high 

sublayer16+j  =  7 

sublayer^  =  6.953  *nm 

xcomppl  =0.15 

pi  =19 

p2  =  38 

dlst  low 

sublayers^.  :=  ? 

sublayerp2  =8.4  -nm 

xcompp2  =  1 

Now  that  all  the  design  considerations  have  been  made,  we  may  start  the  program. 

Incident  angle 

D- 

J5 

*11 

CD 

.  7  upper  “  ^  lower 

ct  :  =  0..  — — - 

SCALE  nm 

^t:  =  Xlower^ctnmSCALE 
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bragg 


^  upper  'K  lower 
nm- 2  SCALE 


X,  850  *nm 

bragg  _ 


nwavez perj0d > ct  :=nAlGaAs(xcomPZperioAt)  ncavZcav>“ :=n  AlGaAs(xcavityzCav  A)  nwave, 

nwave  :  =  Re(  nwave )  o  ncav  :  =  Re(  ncav)  n 


Polarization  selection  made  here.  0  denotes  a  TM  polarization,  1  denotes  a  TE  polarization.  I 
also  give  the  general  expression  for  intrinsic  admittance  of  a  layer. 


adm(n)  :=  Re(n)2-  ( |lm(n)|  )2- 


Pol :  =  1  datag :  =  Pol 


<2/.(  n  y 

Tiadm(n)  =  if!  Pol=0, - ,^adm(n) 

adm  1  adm(n) 


nsubsct :  =  n  AlGaAs(0  00 1  >  O  Set  substrate  to  GaAs  nsubs  : = Re(nsubs)  o 


0in„  =0  Vector  of  all  incident  angles  per  wavelength. 

0,ct 


}  trans  ln 


trails’ n  in 


Function  linking  angle  incident 
from  medium  1  to  angle 
transmitted  into  medium  2. 


0in, 


period’ 


>  trans(nwavez  penod,ct>  nwavez  period_  1>ct,  0inz 


period' 


,ct) 1 


0insubs ,  ct  "  9  trans(nwavelayers  period,  ot.  nsubsct’  0mlayen  period,  ct) » 
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0CaVO,ct  “9miayers  period. ct°  nCaV0,ct  nwavelayers  period>c* 

ecavZcavjCt  trans(ncaVz cav-ct’ ncavz cav-  i,ct.8“vZcw_  1>ct). 


8delay(n  user,  d  user,^  user 


2-n-d 


user 


adm 


n  user) 


'  user 


8z  period. ct  =§delay(nwaveZper.od>ct)  sublayer,  period,  X  J 


Function  describing  phase  delay 
per  finite  sublayer 

Phase  delay  per  sublayer. 


Scayz  cav ,  ct  ■  “  Sdelay  ^ncavz  cav  t  ct,  sublayer  cavityz  *  ^ct 

8topct :  =  8delayfnsubsc[,  20  nm,  L,t 


'•op 

4 


Refnsubs. 


bragg. 


=  58.238  *nm 


This  means  the  layers  add  up  to  it  phase  shifts,  at  the  bragg  wavelength  for  the  periodic  high  -  low 
layers. 


layers  penod 

E 

n=  1 


N 


8  ,  =  1  -  3.737*  10  i 

n,  bragg 


cavity 

8cav 

n=  1 


n, bragg 


=  1.99931  -  0.00277i  *11 


Now  I  will  take  all  of  the  previous  information  and  create  the  Transfer  matrix  described  by  MacLeod  in 
his  book  on  thin  film  filters.  First  I  will  set  up  the  parameters  that  I  use  in  the  calculation. 


Vet  adm(nwaveair,ct)  ^ penod. ct :  =  Vdm(nWaVez  period- * 

ricavz cav,ct  :  =  11  adm(ncavzcav.ct)  Vbs.ct  _T1  adm(nsubsct) 


Tlt0Pct  :  =  11  adm(nsubsct) 


M(8,r|)  :  = 


cos(8) 

j  r|  sin(8) 
siemens 


j  sin(  8)  siemens 
il 

cos(8) 


Characteristic  Matrix  of  layer 


NU=N 


upper 


NU  =6 


NL  :  =  N 


lower 


NL  =  38  LC  :=N 


cavity 


NS  :  =  N  sublayer  NS  =22  LP  :=layers  period  LP  =44  LC“51 
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1  LP  1 

j 

1  LP 

Mupper(a)  :  =  M^8topa>r|topJ- 

n  M(c\j 

i  k=  17  / 

\  k=  1  / 

f  1  ^ 

NU-  1 

1  17  \ 

Mlupper(a)  :  = 

[I  M(Va>\,a) 

\  k  =  LP  1 

a 

ii 

M(8topa>Tttopaj 


1  1  ) 

NL 

1 

23 

Mlower(a)  :  = 

n  m(8m'\.) 

n 

i  k  =  LP  / 

U 

-LP  I 

/  LC 

Mcavity(b)  :  = 

i  I)  M(8CaVk,b,T’CaVk:,b' 

i 

\k=l 

/ 

Now  we  are  at  the  final  page  where  everything  up  to  now  is  put  together.  I  am  using  MacLeod  notation 
throughout. 


‘  mac  ct 


<  >  :  =  Mupper(ct)Mcavity(ct)onoffMlower(ct)- 


'subs ,  ct 

siemens 


calculates  reflectance  of  RCLED 


1 


A1  macV  :  =  Mlower(ct)- 


subs ,  ct 


siemens 


calculates  reflectance  of  lower  mirror 


^macV  =Mlupper(ct). 


'air.ct 

siemens 


'  lower 


mac. 


-t-ct  SCALE 


2,  ct  nm 


calculates  reflectance  of  upper  mirror  from  inside  cavity 


/  \<0>  /  \<0>  /  \<0> 
Bmac  \Amac  )  B1  mac  \A^  mac  )  B2  mac  mac  j 

Cmac:=(AmacT)  siemens  Cl  mac  :=  (a1  macT)  siemens  C2mac:=(A2 


mac 


<l> 


■  siemens 
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Assignations  based  on  MacLeod  notation. 


■'ct.O 


Kct,l 


n  R  —  C 

'•o  ct  mac  t  mac  4 

’  ct  ct 

\  ct  ®  mac  +  ®  mac  , 

’  Ct  Ct, 


^cav5  lct  B*  macc{  ^ 1  macct 
(ricav51  ct-Bl  macct+  macct 


Ticavi  ctB2 mac  ^2 mac 

Ct  Ct 


TR 


4'\ctRe(Wt 


ct,0 


mac  I  ‘  I  ,  ct  B  mac  ,  +  ®  mac 

ct  ct/  \  ct 


TR 


Vet®  mac,  +■  ® 

4TlCaV51,ctRe(Vbs.ct 


"ct,  1 


1lcav5i,ct  ® 1  mac  +■ C1  mac  V  (^^si.ct®1  mac^  +  C1  macct 


TR 


^ct,2 


mac 

ct  ct. 


:’2  (T'cavi  ct  B2mac  +C2 

V  Ct 

4T10,ctRefB 


ct,2 


Ct/  \ 

4TlCaV51,ctReKtf,ct) 


(T1cavijCt'B2mac  +■  C2mac  V  fT'cavi,ct'B2macct‘l"  C2  mac 


AB 


mac .  mac.  Vbs,ct 

ct  ct 


ct,  0 


%,ctB  macc(  +  ^  macctW7V),ct  ®  mac^  ^ 


, .  mac 

ct  ct, 


AB 


4-T)cav  RelBl  mac  Cl  mac  Vta.et 

*  \ _ Ct _ Ct _ 


Ct,  l 


(T|cav5i  ct'Bl  mac  .  +‘  C1 


AB 


mac^1-  C1  macj  ^caV5i,ct'B1  mac£t  +  C1  macct 
4  r)cav51  ct-Re|B2  mac^  C2  macct  ~  Vr.ct 


ct,  2 


mac 


ct, 


|Tlcavi ,  ct  B2  macct  +  C2  macj '  ^cavi ,  ct  ®2  macct  +  C2 

Complex  reflectivity  of  field  derived  by  matrix  method 

3hase  change  upon  reflection  by  DBR  stack. 

q  =0..  2  Vq  ^  (angle(Re(Pct,q)>Im(Pct,q)))  RCt,q^l 


P 

q  ct,q 


R  :=4 

ct,q-t-4  ct,  q 


Rct, o  Amac2>ct  j  :=0..7 


R.  :  =  data. 

J.7  j 


R  .  „  =TR 
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Intensity  reflectance  of  DBR  Stack 


Useful  to  have  all  user  input  parameters  right  here  beside  the  graph 

Mower"750  11111  7  upper  ~050  nm  Range  of  wavelengths  used  by  program  Po1  1 

SCALE  =  1  kop  =  850nm 

0=5  deg  N  upper = 0  ^  lower onoff=l 

J  .  .  .  WRITEPRN(Ref405)  :  =  Ro 

x  lowjndex  2 1  0  Composition  of  lowest  index  desired 


x  highmdex"015  Composition  of  highest  index  desired 


The  next  step  is  to  use  the  information  previously  calculated  to  determine  the 
electric  field  intensity  at  each  layer  in  the  structure.  The  first  thing  to  do  is  create 
the  entire  RCLED  structure  by  stacking  the  matrices  with  index  and  phase 
information. 

N  sublayer  =  ^  laYers  period  ~~  ^4  N  iower  -  38  N  Upper  -  6 
^  cavity  ^  sublayer  -  75  ^  cavity  -  5  ^ 

1 1  :  =  layers  perj0Cj  N  Upper N  cavj(y+ layers  perj0(j  N  jower  1 t  -  1.987  10 

1  ]  -  layers  penocj  N  iower  j  _  j  572*  jq3 
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1  u  =  264 


1  u  layers  penod'  N  Upper 
*  c  '  ”  ^  cavity  ^  c  ~  3  ^ 

rows(nwave)  =45  cols(nwave)  =201  rows(ncav)  =  52  cols(ncav)  =201 

z  :  =  0..  3 

mak  1  1.1  u  n0  z  1  choice  :  =  ( bragg -  25  bragg  bragg- 8  bragg-40)T 

bragg  =  100 


nmakl ,  z  =  if(m0d  (makl  ’ layerS  period)  *°  ■  nWaVelayers  period,  choicez>  nwavemod  (makl ,  layers  period) ,  choice.,) 

^makl  :  -  if(mo<l  (mak  1 , layers  period)“0’®inlayersperiod,choice>6inmod('mak1, layers  period), Choice)  0 

SUbmakl  =  if(mod(m£d<1  >  layers  period)*^  >  SUblayer|ayers  period’  SU^^ayermod (makl, layers  period)) 

mak2  -  1  -t- 1  u--  N  cavjty  -l- 1  u 

92mak2  -9caVrnod(mak2,l  u)  .choice- 2°  nmak2,z  ~ nCaVmod(mak2, 1  u) ,choicez 


submak2  :  sublayer  cavity^d^u,!  u) 


62o:  =  0 


mak3  :  =  1  +- 1  c  -t- 1  u~  1 1  rows(nwave)  =45  i  :  =  1 ..  layers  period 


nwave2. 


1,Z 


nwave 


<(  .  ^ 

choice 


sub2.  :  =  sublayerlaye„per.od  +  1_i 


J  layers  period  H-1-1 


. .  -  ,  <choice  — 2> 

0m2.  :=  i0in 


/layers  period  d-t-i0 

nmak3,z  =  ^(m°d  (mab3  “  1  c  >  layers  period) *°  ’ nwave2layers  period. z’ nwave2mod(mak3  -  1  c.  layers  penod)  ri 


92mak3  :  =  (m°d (mak3  ‘  1C’ layerS  period) *° ’ 9m2lavers  penod=  9in2mod (mak3 -  1  c. layers  peri0d) )  ■ 


SUbmak3  :  =  lf(m0d  (mak3  -  1  c  > layerS  period)  *°  ’  Sub2layers  period  >  Sub2mod  (mak3  - 1  c.  layers  period)  ) 


rows(n)  =  1.988- 103  1 1  =  1.987-103 
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Yeh  method  from  "Optical  Waves  in  layered  Media' 


D(n)  :  =  ifl  Pol*0 , 


1  /  adm(n)  adm(n)\  /  1  1  \ 


n  \  1  -1  /  \adm(n)  -adm(n) 

i  k  v  d 


P  kx»d  - 


\  0  e 


0  |  2-7! 

k(n,X)  :=—  (adm(n)) 


i  ;  =  0..  1  countl  :  =  1 1 -  2..  0 
/O.OOl  \  volt 


AMPi>lt_!  :  = 


AMP1 


0  .  m 

1  1 


Mt-l  " 


0.00 1  \  volt 


0  /.  m 

1 1 


use  :  =  0 


AMP<COun  “P  k[ "countl , use ’ k (choiceuse) ] ’ SUbcountl  D  (“countl,  use)  D (“countl •+- 1, use)  AMP 

AMP<0>  :  =  D  (n  W-Dfn,  ')  AMP<1> 

1  0,u scj  V 


,<countl  -f- 1> 


use ;  =  2 


AMP1<C0Unl  ' “ P( k(  “count  1 ,  use ’ k (choiccUSe)  ’^countl  D  (“countl, use)  D  (“countl -(-1,  use)  ’ 


,<countl  -\- 1> 


AMP1<0>  :=D(n  W-Dfn.  )-AMPl<1> 

1  0,u sey  V  l,usc/ 


choice1  =(  75  100  92  60  ) 


z total  "  1  " 1 1  “  1 


P"z  total  AMP0.ztotal  +  AMPl»z  total 
P'  1 z  total  =  AMP 1 0,  z  total  AMP 1 1 ,  z  total 


Calculated  E-Fields 


^  total  =R^ 


lt-l 


sub 


cl_z  total 
1 


Use  to  reference  each  sublayer  to  thickness 


I“te“sityz  to^ , o  ^  "^,0  ^  (e^.  E,  J  max( Intensity)  =5.829- 10"5  --Watt 


m 
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m 


Mensit^  total,  i  :=2IW3-fl'lE1WE1*totJ  max  (intensity''1  > )  =3.385- 10^4 
norm  f : = max(  Intensity ) 

use  =2 


Re 

/ - >\ 

Intensity 

nornif  j 

norm  j  :  =  max(n) 


V2z total-0  :  "z total- ^  3 


V2zt„t,,l  "z 


total’ 1  ~  total 


WRITEPRN(ELEC5 )  : = OUT  o 


Appendix  C.  Quantum  Well  Energy  Level  and  Gain  Calculations 

This  document  calculates  the  energy  levels  of  rectangular  quantum  wells  in  bulk  material.  The 
quantum  wells  are  finite  and  have  the  confinement  energy  defined  as  V0. 


kb  =  1.380658- 10"  23-i^ 
D  K 


eV 


=joule  1.60217733- 10' 


.-19 


nm=mT0"9  Hm^m-10 6 


hbar : 


.6.6260755-10' 
2-n 


.-34 


■sec  h:  =  hbar-2-7i  X(E)  :-h— 
joule  E 


=8.854187817- 10" 12:^ 


m 


q  =1.60217733- 10" 19  coul  me  =  9. 1093897-  10"31-kg 


c=299792458- 


m 


sec 


The  preamble  defines  all  the  terms  that  will  be  used  in  the  document.  User  entered  parameters 
follow: 


m  re(x)  =0.067  +■  0.083  X  Relative  mass  of  electron  as  function  of  composition  x. 

m  rlhW  =  0  087  +■  0  063  x  Relative  mass  of  light  hole  as  function  of  composition  x. 

m  rhh(x)  =0.62  +  0. 14  x  Relative  mass  of  heavy  hole  as  function  of  composition  x. 

x  well  =0  0  Composition  factor  of  Quantum  well.  Used  to  calculate  conduction  and  valence 
band  energies 

x  cav  :  =  0. 1 5  Composition  factor  of  cavity 

L  =8  nm  Thickness  of  Quantum  well  in  nanometers 


me  well  m  re  (x  well) ' m  e 

Mass  of  electron  inside  quantum  well 

me  well  =  0067 ’me 

me  cav  _  m  re(x  cav)  m  e 

Mass  of  electron  outside  quantum  well 

me  cav  =  0  079  *me 

m^1  well  • -  m  rhh(x  well)  m  e 

Mass  of  heavy  hole  inside  quantum  well 

“Mi  well"  062 *me 

m^1cav  ~mrhh(xcav)  me 

Mass  of  heavy  hole  outside  quantum  well 

mhh  cav  =  0.641  *m  e 

well  m  rlh(x  well)  m  e 

Mass  of  light  hole  inside  quantum  well 

mih  wejj  =  0.087  *m  £ 

m^1  cav  “  m  rlh(x  cav)  m  e 

Mass  of  light  hole  outside  quantum  well 

mill  cav  =  0.096 ’m  e 

For  A^Ga^xAs 
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Eg(x)  :  =  if(x<0.45, 1 .424  +  1 ,247-x,  1 .9  +-  0. 125-x  +-  0. 143-x2)  eV 

ECl(x)  :=Eg(0)  h- 0.6  (Eg(x)  -  Eg(0))  EVl(x)  :=-0.4(Eg(x)  -  Eg(0)) 

Used  method  of  Bour  to  come  up  with  following  graph 

xdemo  °5  -  1  Range  used  to  create  a  sample  graph 

The  composition  range  where  the  conduction  band  energy  (referenced  to  the  0  composition  zero 
point)  is  rising  is  where  the  semiconductor  is  a  direct  bandgap  type,  k  conservation  rules  apply  here. 


I  will  now  define  the  characteristics  of  the  quantum  well  using  the  input  parameters  above.  From  the 
compositional  fractions  of  the  well  and  cavity  I  will  create  an  energy  diagram  for  the  cavity.  I  am 
using  A^Ga^s  as  the  cavity  and  well  material. 


EC  cav(xl,x2,d) 


t,ECl(xl),ECl(x2)J  EV 


cav(xl,x2,d) 


<^,EVl(xl),EVl(x2)J 


d  :  =  -  10nm,-9.91nm..  10- nm 


EC  cav  (x  wc]| ,  x  cav,  d) 


eV 


6i  EV  cav  (x  Wclbx  cav  ’  4) 

5 


Energy  Diagram  for  Cavity 


Distance  across  Cavity  (nm) 
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V0  is  defined  as  the  difference  between  the  energies  of  the  conduction  bands  between  the  quantum 
well  and  the  cavity.  V0  is  zero  referenced  to  the  top  of  the  conduction  band  of  the  well  material. 
There  are  three  V0s  defined.  One  for  electron  energy  levels  in  the  conduction  band  of  the  quantum 
well,  and  two  for  the  Ih  and  hh  energy  levels  in  the  valence  band. 

V  %  :  =  ECl(x cav)  -  ECl(xwdl)  ECl(xcav)  =  1.536-eV  ECl(xwell)  =  1.424 -eV 
V0i:  =  EVl(xwell)-EVl(xcav) 

Vo2  :  =  EV1(X  well)"  EV1(X  cav)  V  QT  =  (  0.1 12  0.075  0.075  )*eV 

k0  is  a  constant  used  in  the  calculation  process.  It  simplifies  the  bookkeeping  by  wrapping  up  a  lot 
of  other  numbers  into  one.  It  has  units  of  inverse  meters. 


k0T  =  (  0.408  1.085  0.393  )*ran  1 

Now  I  use  the  secant  method  to  solve  the  transcendental  equation  posed  by  the  electron  energy 
levels  inside  the  quantum  well.  The  equation  is  highlighted  and  explained  in  the  book  by  Weisbuch. 

yl(x)  :=x  tan(x)  y2(x)  :=xcot(x)  y3(x,r)  :=^r2- x2 

f(x,r)  :  =  yl(x)  -  y3(x,r)  g(x,r)  :=y2(x)-i-y3(x,r) 

Because  the  equation  has  different  forms  depending  upon  whether  you  are  looking  for  an  even  or  odd 
quantum  level,  I  must  ensure  that  Mathcad  knows  which  condition  to  apply,  r  is  a  tidy  constant  to 
wrap  things  up  a  little  and  help  the  form  of  the  equation,  r  is  dimensionless. 

L  me  well  ,  L  “^well  ,  L  "^well 

V=  V=  °i  2  mhh  cav  ^  "  °2  2  mlhcav 

rT=(  1.376  4.199  1.416  ) 

From  r  we  get  the  number  of  allowed  energy  levels  in  the  quantum  well.  Once  we  know  how  many  we 
have  to  solve  for,  it  becomes  easier  to  find  them. 
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Now  I  apply  the  secant  method  to  solve  the  equation.  I  know  that  there  are  sol  roots  so  I  set  the  root 
counter  n  from  0  to  sol  -1 .  The  secant  method  requires  two  initial  guesses.  These  may  be  varied  as 
necessary. 


n  :  =  0..  sol0  -  1 

q  :  =  0..  solj  -  1 
p  :  =  0.. sol2-  1 


Initial  Guess  1 


u  71 

xc  :  =  —  +  n— 

n,o  22  2 


n  n 

xv  :  =  — -t-q— 

q,°  2  2  2 


.  71  71 

xvl  :=  —  +  p— 
P,°  2  2  2 


Initial  Guess  2 


n  7t 

xc  :  =  —  +  n— 

n’1  2.7  2 


71  71 

xv  :  =  — +  q— 
q>1  2.7  2 


,  71  71 

xvl  :  =  —  +  p - 
P>1  2.7  2 


N  denotes  how  many  iterations  to  use  to  converge  to  a  solution.  My  initial  guesses  already  comprise 
the  first  two  solution  "slots". 


N  - 12  i  :  =  2..N 

The  solution  variables  change  depending  upon  whether  an  odd  or  even  solution  is  being  sought.  The , 
odd  or  even  solutions  to  the  wave  function  represent  cos  or  sin  wave  functions  across  the  well.  (Ref 
Weisbuch) 


f(xl,x2,r,n)  :  = 


mod(n,2)ss0,xl  -  r> (xl  -  x2)  xi  _  |(^,rHxl^2) 

f(xl,r)-f(x2,r)  g(xl,r)  -  g(x2,r)/ 

T 


<N>  ■ 


=  0.878 


xc  .  :=f(xc  ,,xc  .  Electron  Solution  xc 

n,i  \  n,i-l  n,i-2  0  J 

xv  :=f(xv  .  ,,xv  .  _,r,,q')  Heavy  hole  solution  xv^  =(  1.265  2.503  3.656  ) 

q,i  \  q,i—  l  q,i-2  1  V 


xvl  .  :  =  f^xvl 


p,i-  l’ 


,xvl 


p,i  — 2  2 


.r,,p 


Light  hole  solution  xvl 


<N> 


=  0.891 


The  following  graph  plots  the  calculated  points  at  the  intersections  of  the  transcendental  equations. 
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xc  :  =  0.01, 0.02..  max(r) 


The  answer  derived  before  is  of  the  form  x  =  kU2.  To  get  the 


kc  :  =  xc 


<N>[  2 


kcT  =0.22  *11111  1 


kv  :=  Ixv^l--  •  kvT  =  (  0.316  0.626  0.914  >nm_1 
L 


allowed  k  from  this  solution  I  must  use: 

Electron  allowed  wavenumbers 

Heavy  hole  allowed  wavenumbers 


kvl  :=  xv  1 


<N> 


kvl T  =  0.223 -nm'1 


Light  hole  allowed  wavenumbers 


The  relationship  between  k  and  Energy  for  particles  is  not  linear  but  is  related  by  the  particle's 
momentum  p.  The  relationship  is  given  by  (for  the  quantum  well). 


^  well  (k  ’  ^  o  ’ m  well)  ■ 


hbar2k2 ' 
^'m  well/ 
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Now  I  put  the  k  solutions  back  into  the  energy-k  equation  and  solve  for  the  allowed  energies.  I  will 
then  reference  them  back  to  the  original  zero  reference  defined  by  Lott. 


EC  aiiown  :  _  EK  well  ^cn  >  V  oQ  > me  well j  EC  auow  0.027  *e V 

EV  =EKwell(kv,,V0i,mhhwell)  EV  ,„0WT  =(  0.006  0.024  0.051  >eV 
EV1  allow  =EK„||(kvlp  V0j.mlhwd|j  EVI  allowT  =  0.022'eV 

Now  I  reference  these  back  to  the  original  energy  levels 


EC  orig  EC1  (x  well)  +■  EC  auow 

ECorigT  =  l.451-eV 

EV  orig  "  E V1  (x  wen)  -  EV  anow 

EvOrigT=(-0.°06  -0.024  -0.051  )*eV 

EVI  orlg  _EVl^xwej])  -  EVI  auow 

EVI  origT= -0-022 -eV 

EVTorig:  =  Sttdc(EVori8,EV,ori4 
on8  \  eV  eV  )  eV 

EVT  origT  =(-0.006  -0.024  -0.051  -0.022  )*eV 

Now  I  have  to  define  the  transition  energies  that  are  possible.  This  basically  means  finding  the 
difference  in  energy  level  between  every  electron  level  and  heavy-hole,  light  hole  level. 

s  :  =  0..  ^solj  +  sol2)  -  1 
E  trans  ~  EC  or,g  -  E  VT  ong 

n,s  n 

E  =  (  1458  1475  1503  1473  >eV 

- > 

E>,(E)  :  =  — -  l1111  ^-trails  -E^(E trans) 

E  ev  v  ' 

51  trans  =  (  850  731  840.401  825.151  841.739  >nm  Transition  Wavelengths 

Now  I  will  apply  what  I  have  already  done  to  the  Chinn  model  to  calculate  the  gain  and 
spontaneous  emission  for  a  cavity  with  no  mirrors. 

Mas(  |M|  )2  From  Chinn  and  Lott  notation,  this  is  transition  matrix 

E  g  :=ECl(xwell)  -  EVI (xweU)  E  g  =  1-424 -eV 
A0(x)  :=(0.34-  0.04- x)  eV  Ao(0)  =0.34'eV 
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M  =  1.893  •eV  me  Transition  Element  Matrix 


reduced_mass  fm  j ,  m  2)  :  ~  l"  1  +  m  2~  lj 

m  electron  :  =  mewell  m  electro/ =  0067 ‘m  e  Reduced  mass  electron 

n 

mjjjjQjg  :=mhh  well  m  hhole1  =  (  062  0  62  0  62  )*me  Reduced  mass  heavy  hole 
q 

mjjjQjg  :  =  mlhwell  m  moleT  =  0.087  *m  e  Reduced  mass  light  hole 

P 


M  :  = 


me2Egi 

VE  g"*"  ^o(xwell)) 

1 2' me  well'  1 

^o(xwell/ 

E  o-i- 2-  v  y 

l  g  3  / 

mhole  :  =  stack 


m  hhole  m  lhole 


m , 


m. 


■m , 


m 


hole 


mtransl  “re^uce<^-massfm  electron  > m  hhole  ) 

n,q  \  n  q/ 

mtrans2  =  reduced_mass/m  eiectron  ,m  lhole 


m  trans  augment  l  m  ^  j ,  m  trans2 


=  (  0.62  0.62  0.62  0.087  )*m  e 
m  trans  1  =  (  0  06  0  06  0  06  )'me 


mtrans2=0038‘me 


Now  that  I  have  defined  all  of  my  transition  energies  and  reduced  interaction  masses,  I  will  define  the 
selection  rules.  I  am  initially  enforcing  rigorous  k-selection.  Only  transitions  between  like  quantum 
states  are  allowed.  That  is,  1  -1,2-2,  etc. 

Cln  :=if(n=q,  1,0)  Cl  =(  1  0  0)  C  matrix  for  electron  -  heavy  hole  transitions 

C2  =  iff  n=p  ,0,0)  C2  =0  C  matrix  for  electron  -  light  hole  transitions,  however  we  are  only 
n,p  considering  electron-heavy  hole  interactions,  therefore  this 

parameter  is  always  zero 

C  :=augment(Cl,C2)  C=(l  0  0  0  )  Complete  C  selection  matrix 

El  \  Anisotrophy  factor  for  the  transitions,  considering  only  TE 
*  e2  I  radiation  for  the  electron  -  heavy  hole  transition 


A(E1  ,E2)  :=— • 
4 


fcv(E»Ef.T) 


1 


lE-Ef\ 

expj 

+ 

\  kb' 1  1 

1  1 

Fermi  population  factor 
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Fermi  population  factor 


fcv(E,Ef,T  := 


/E-Ef\ 

exp 

ikbTi 

+- 

1 

H(y)  :=if(y>0  eV,  1 ,0)  Heaviside  Step  Function  used  to  define  allowed  energy  transitions 
First  I  must  find  the  quasi  Fermi  levels  for  the  concentrations  and  temperature  specified 


T  :  =  300  K  E  f  :=  1.4605-eV 


carrier  :  =  (1  1.5  2.0  2.5  3.0  3.5  )T  1018 


EC  orig  "  1  451  -eV 

Assuming  Charge  Neutrality  q  :=0.,5 


CONC  eiectj-on  (E  g.  j  log 


kb'T 

hbar2Jt-L 


Sm  electron  k1 

n 


1  +-  exp 


EC  orign  E  fc 

kbT 


_  _ 

- 

3 

■cm 

CONChole(EfV)  ;  =  log 


kbT  ■ 

hbar2-  n-  L ' 


“hole  “ 


1  -i-  exp 


Efr-EVT 


ong. 


kbT 


■  i 

■ 

3 

•cm 

conc electron  :  =  l°g(  carrier)  conchole  :=1°g(' carrier) 


i 

guess £c  :=  1.51-eV  guess^,  :  =  0.001  eV  AEg  :  =  -310'8  eVcarrier3 


E  fc  -root/ CONC  electron  (g1*83  fc)  "  conc  electron  ’^ess fc 

q  \  '  '  1 


Efv  -root (CONC hole (guess A  -  conchole  ,guesSfV' 
q  \  q 


TOL  =  10' 


12 


ml 


E  e(El  ,ml  ,m2)  :  =  EC1  (xwdl)  ^  (El  -  E  g) 

Eh(El)ml,m2):=EVl(xwell)-(El-Eg)m2 

Finally,  put  it  all  together  for  the  gain  spectrum. 

qC2M 


m2  Conduction  band  energies 

ml  Valence  band  energies 


con  :  = 


, ,  „  _  T  Constant  defined  for  calculations  _ 0  r2 

£0me2'Chbar-3.6L  con  =  1.06*10  *m*sec 
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Etransl  ■  ^trans  -'"AEg 
q  0,0  Dq 

EtranslT=(  1428  1423  142  1417  14,4  1412  >eV 

Fermi  Inversion  Factor  given  by 


fg(E,z)  fcv^E e(E’m transQ  0,m  electron^ ,E  fy  4 

Probability  of  Emission  given  by 

fe(E,z)  -fcv^Ee^E,mtrans0  0,melectron0j,Ef<y  4j' 

Gain  spectrum  given  by: 


^cv^Eh^E,mtrans0  0’mhhole0j’Efy4j 

1  “  ^cv^Eh^E>mtrans0  0,mhhole0j,Efvz’4Jj 


g(E,z)  ■con-mtrans^  ^  A^E,Ejransj^  fg(E,z)-H^E  E^^jJ 

Now  convolve  the  calculated  gain  with  the  Lorentzian  to  account  for  intra-band  scattering 


i  :  =  0..  1023  Etest  :=  1-eVq- i-.002  eV 
i 


j  =0  1023  Xc  :=  10" 13  sec  Numberfrom  Corzine  (Zory:  Quantum  Well  Lasers)  accounts 

for  Lorentzian  broadening 

gltest^q  : =  g^E  test.’  qj 


k  =0..  511 


hbar 


lor(E)  :  =  — 


71  L  max(E  test)  +  ^'eV\ 


Lorentzian  used  for  intra-band 
scattering 


Lor.  :  =  lor(Etest 


•eV  Lori  :  =  Lor 


L0rk:  =  L°rl512+k 


Lor, 


512-t-k 


:Lorl, 


Now  normalize  convolved  spectra  to  original  spectra  and  print 


Gtest<q>  :=iffl((ffl(Lor)-ffi(gltest<q>) ])  Gltest<q>  :=  Gtest<q> ■  ^ 

iGtesH 
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Again  convolve  with  Lorentzian  and  normalize 


;i,q  : ■  R(E  test.’ 3 j  G2test<q>  := ifft((fift(Lor)  fft(Rs<q>)  j)  Rltest<q>  := G2test<<1> 


ZRs<q>  ZG2test<q>  lRltest<q> 


5.291-  1054  kg~1m“5sec 

8.234-  lO^-kg^-m^-sec 

5.291-  lO^-kg’^m^-sec 

1.086-  1055-kg‘1-m'5-sec 

1.69-  1056-kg‘1-m“5-sec 

1.086- 1055  kg‘1-m‘5-sec 

1.785-  1055-kg"1-m'5-sec 

2.778- 1056-kg_1-m'5-sec 

1.785- 1055  kg-1  m“5-sec 

2.556-  lO^-kg*1  nT5-sec 

3.977-  1056-kg'1-m"s-sec 

2.556-  1055-kg"1-m"5-sec 

3.414  1055-kg”I-m"5-sec 

5.313-  1056-kg‘1-m’5-sec 

3.414-  lO^-kg'^m^-sec 

4.35-  1055-kg‘1-m'5-sec 

6.77-  1056-kg"  1  m’5-sec 

4.35-  1055-kg'1-m"5-sec 

carrier1  =  (  1  1.5  2  2.5  3  3.5  )*1018 


lRs<q> 

lG2test<q> 


Next  few  lines  output  data  to  file 


result.  „ 
1,0 


nm 


3 

result.  ,  :  =  Rltest.  •  eV-cm  -psec 

1,1  1,1 

result.  _  :=Rltest.  ■  eV-cm3 -psec 

1,2  i,2 

3 

result.  :  =  R 1  test.  3- eV-cm  psec 

result.  .  :  =  Rltest.  •  eV-cm3  psec 

1.4  i,4 

result.  , :  =  R 1  test.  ■  eV-cm3  psec 

1.5  i,S 


carrier  =  2.5x1 018 
carrier  =  3.0x1 018 
carrier=  3.5x1 018 
carrier  =  3.0x1 018 
carriers  4.5x1 018 


Put  write  equation  on  toggle,  so  it  only  writes  when  you  want  it  to. 

WRITEPRN(rspon)  :  =  Re( result)  o 
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